Wavefront sensing system employing active updating of reference positions and subaperture locations on wavefront sensor

ABSTRACT

A free-space adaptive optical laser communication system having signal transmission and reception channels at all terminals in the communication system, wherein wavefront sensing and wavefront correction mechanisms are employed along signal transmission and reception channels of all terminals in the communication system (i.e. adaptive optics) to improve the condition of the laser beam at the receiver (i.e. reduce the size of the spot the detector plane). Speckle-to-receiver-aperture tracking mechanisms are employed in the transmission channel of the communication system and laser beam speckle tracking mechanism in the reception channels thereof, so as to achieve a first level of optical signal intensity stabilization at signal detector of each receiving channel. Speckle-to-fiber/detector locking mechanisms are also employed in signal receiving channels of all terminals in the communication system so as to achieve a second level of optical signal intensity stabilization at signal detector of each receiving channel.

RELATED CASES

The present application is a continuation-in-part (CIP) of co-pending U.S. application Ser. No. 10/777,476 filed Feb. 12, 2004; said patent application is assigned to and commonly owned by Metrologic Instruments, Inc. of Blackwood, NJ, and is incorporated herein by reference in its entirety.

BACKGROUND OF THE INVENTION

1. Field of Invention

The present invention relates to an improved bi-directional laser communications link between two optical transceiver platforms (e.g. an airborne platform and a satellite) having a way of and a means for dynamically stabilizing variations or fluctuations in the intensity of detected optical carrier signals (i.e. suppressing fade events) that are caused by atmospheric effects such as turbulence, scintillation etc. which tends to degrade link performance. Such fade compensation is particularly important for free-space optical laser communication systems requiring the transmission over minimum laser beam power levels over fiber-to-fiber free space optical communication links.

2. Brief Description of the State of the Art

There is great interest in super broad-band free-space optical (FSO) laser communication systems as they are capable of securely communicating information at high data rates in point-to-point and multi-point-to multi-point communication networks.

In military applications, such laser communication systems and networks offer a level of superiority and security over radio-frequency (RF) based communication system which have relatively limited band-widths, and thus data transfer rates, as well as being susceptible to RF-based jamming techniques intended to interfere and disrupt the performance of such systems. In commercial applications, such laser communication systems can be rapidly installed in point-to-point and multi-point-to-multi-point configurations (using buildings and towers as support structures for such laser communication platforms) at a significantly reduced expense in comparison with micro-wave-based satellite communication systems. Examples of such FSO laser communications systems are disclosed in U.S. Pat. Nos.: 6,657,783; 6,643,467; 6,348,986; 6,347,001; 6,314,163; 6,285,481; 6,286,944; 6,181,450; 6,151,340; 6,122,084; 5,844,705; 5,786,923; 5,786,923; 5,754,323; 5,710,652; 5,606,444; and WIPO Publication No. WO 03/003618, each said prior art reference being incorporated herein by reference.

However, FSO laser communication systems are not without challenges and problems and challenges.

In particular, the free-space optical communication links in FSO laser communications are effected by atmospheric conditions such as turbulence and the like which works to cause aberrations in the spatial phase of the wavefront of the modulated carrier laser beams as such carrier laser beams are transmitted between terminals in such communication systems. Depending on the distances between the laser communication terminals, such spatial phase aberrations evolve into spatial intensity aberrations in the laser beam received at the entrance pupil of the receiver module of such terminals, degrading the bit error rate (BER) of such communication systems.

Studies have been conducted by Adaptive Optics Associates, Inc. to quantify these wavefront errors over the range of possible link geometries. The uplink and downlink have been considered separately because the atmosphere affects the propagation differently depending on its direction. The density of the atmosphere declines rapidly with altitude and thus most of the wavefront errors are generated close to the aircraft. Integration of the analytical model of the atmosphere proposed by Fried (JOSA, 56,1380, 1966) shows that, for the altitudes of interest in this study, the effective distance of the turbulence from the aircraft is about 3 km.

For the light propagating up to the aircraft the wavefront errors picked up close to the aircraft evolve along the long path to the satellite so that the beam that reaches the satellite has both phase and intensity errors. This behavior is described by the eikonal equation (Born and Wolff, Principles of Optics, p. 112), which shows that the change in intensity between to planes perpendicular to the propagation direction is proportional to the product of the curvature of the wavefront and the propagation distance. This effect is the source of scintillation in the beam.

For light coming down to the aircraft, the situation is very different. The wavefront has only 3 km to evolve before reaching the receiver and thus is almost purely a phase error without any scintillation. Therefore, the airborne to spaceborne direction as being the most problematic direction in this link. This fundamental difference leads to very different requirements for the satellite and aircraft transceiver compensation systems.

Modeling of Atmospheric Turbulence

To assess the impact of atmospheric turbulence on the communications link, a model of the atmosphere has been developed. This model is based upon one of the standard models (e.g. Clear 1 model developed for the ABL program) supplemented with field data from the ABL-EX experiment, data from astronomical observatories (Masciadri, et al. Astron. Astrophys. Suppl. Ser. 137,185-202, 1999.), and the AFRL 3D modeling project. Of particular importance in this study is the modeling of the turbulent layer associated with the tropopause, determination of the range of altitude for the tropopause, and modeling of the effects of a jet stream. This is because, unlike an observatory that can be sited at a location known to have very low turbulence, the communications link must operate even if the aircraft is below the tropopause and a strong jet stream is traversed at some zenith angles.

Using the atmospheric, turbulence profiles along various link paths have been determined. The important parameters are the altitude of the aircraft, the velocity of the aircraft, and the altitude and orbital path of the satellite. The aircraft altitude determines the starting point of the link path. The aircraft velocity sets the rate at which the beam translates through the atmosphere. The satellite altitude and path determine the apparent slew rate of the beam. FIG. 1 describes these platform parameters. Once a set of link paths is defined from these parameters, the atmospheric model has been used to create a series of phase screens at different altitudes. These are inputs to the AOA Wave Optic Propagation Code that is used to calculate the effects on the laser beam. Such modeling has been performed for both the uplink and downlink directions. One additional input to the Wave Optic Propagation Code is the aero-optical effect of the flow field around the window of the transceiver on the aircraft. This is modeled as a time varying phase screen. The results of this simulation will be a time sequence of phase and intensity maps for the entrance apertures of the two receivers.

In FIG. 1B a typical phase map calculated for a horizontal communications link with strong turbulence. FIG. 1C shows a time sequence of detected laser beam spot intensity distributions calculated using a simple model atmosphere for propagation to a satellite at 1000 km altitude from an aircraft at 35,000 feet. The bar in each frame represents 10 m at the satellite. The second through fourth frames in this time sequence show a dark (rather than bright) speckle of several meters extent. Such an event would lead to a deep fade of the link signal for any receiver aperture of less than a few meters diameter.

The atmospheric model described above is used to examine other aspects of the laser communication link compensation problem. For example, because of the long path to the satellite, the light travel time to the satellite becomes significant. The result is that the transmitted beam must be pointed ahead of the apparent location of the satellite, just as a hunter must point his gun ahead of a flying duck to successfully hit it. This point-ahead angle is typically tens of microradians. Difficulties arise when this point-ahead angle becomes larger than the isoplanatic angle of the atmospheric turbulence. This is the angle over which the aberration is correlated. If the point-ahead angle is larger than the isoplanatic angle, the wavefront error measured from a source co-located with the satellite does not correctly represent the aberrations that will be experienced along the path of the transmitted beam. Also, as wavefront tilt is the strongest atmospheric aberration, this aniosplanatic effect can lead to serious pointing errors in the transmitted laser beam if the employed laser beam tracking algorithm uses only the satellite apparent position to estimate the error. A great body of research exists on this effect, both in the astronomical literature and from HEL DEW projects. If tilt anisoplanatism appears to be a significant problem for a given application, then more advanced tracking algorithms need to be drawn from this prior research and tested in the simulation.

Using the above-described techniques, some initial modeling has been performed to verify that atmospheric turbulence is a significant issue for aircraft to space optical communications links. A model atmosphere was constructed based on the analytical model of Fried (JOSA, 56,1380, 1966). The Fried model was modified to include the observed “bump” in turbulence in the region of the tropopause and the strongly layered structure of turbulence. FIG. ID shows a typical C_(n) ² profile generated by this model. This model was used to estimate the integrated turbulence along a vertical path from aircraft at various altitudes. The results for one realization of turbulence are shown in FIG. 1E. From these integrated turbulence estimates pertinent parameters for the design of a compensation system were derived. These included the atmospheric coherence length (Fried's r_(o) parameter) and the coherence angle, known as the isoplanatic angle. These calculations were made using a wavelength of 1550 nm.

FIGS. 1F1 and 1F2 show some typical results derived from modeling this phenomenon. In particular, FIG. 1F1 shows that for altitudes below the tropopause (about 45,000′ in this case) the coherence length is comparable to, or smaller than potential transceiver aperture sizes. This means that significant atmospheric aberrations will occur at these altitudes. FIG. 1F2 shows that at those same altitudes, the isoplanatic angle is in the range of tens of microradians. Since this is comparable to the point-ahead angle required for a link to a LEO satellite, the issue of tilt anisoplanatism must be treated carefully in any particular application.

Wave-Optic Propagation Simulations

To further assess the impact of this level of turbulence on a free space optical communication link, some simple wave-optic propagation simulations were performed. Using the integrated C_(n) ² value, a single phase screen was generated for an aircraft altitude of 35,000 feet. FIG. 1G shows a pseudo-color representation of such a phase screen. The full extent of the screen is 1m×1m. Subsections of this phase screen were extracted along a path to simulate the motion of the aircraft past the turbulence. In this simple model, the turbulence itself is frozen in time. The extracted phase screen was then applied to a uniformly illuminated circular pupil of 150 mm diameter and the resulting wavefront propagated to a range of 1000 km.

While the use of a single phase screen does not model the real three-dimensional distribution of the turbulence, because the atmospheric density falls off fairly rapidly with height, it is adequate for obtaining a rough idea of the effects of turbulence on the beam pattern at the satellite. FIG. 1H shows two examples of the beam profile at the satellite compared with the diffraction limited profile (left image). Each image covers about an area 100 m by 100 m. In the case of the center beam profile, little power is lost due to the aberration. However, the laser beam profile on the right shows that significant fade events can occur as the receiver aperture could easily lie within one of the dark (speckle) regions of the laser beam spot image.

A sequence of 576 beam profiles was generated during wave optic propagation simulation. For a platform velocity of 200 m/s this sequence represents about 0.1 s of flight time. The temporal variation of the received power was calculated for this sequence and is shown in FIG. 1I. There is one fade event of about 10 db in this data set. Thus, in any link compensation system, it will be necessary to determine how well it suppresses these occasional fade events in addition to the expected improvement in link power due to the improved beam pattern.

The effect on laser beam pointing has also been examined. For each profile, the location of the centroid of the light was calculated. The rms variation in pointing direction was found to be 12.6 μrad with a peak deviation of 28.7 μrad. This rms variation in pointing is slightly larger than the diffraction limited beam radius and thus must be corrected to maintain an adequate communication link signal.

To determine the parameters for a compensation system, the spatial frequency content of the phase aberrations must be examined. Each of the extracted phase screens is decomposed into Zernike polynomials and the temporal intensity fluctuations of the Zernike coefficients determined. FIG. 1J displays the rms variation in each of the lowest 15 Zernike terms. Notably, the bulk of the aberration is contained in the first three terms, the two tilt terms and defocus. That is typical of all atmospheric aberrations. Correction of just these three terms, in transmission and reception channels of the system, would produce a laser beam that, in a statistical sense, would have a power density close to the diffraction limit. For typical imaging or directed energy applications, this is the figure of merit used to assess the degree of compensation to be achieved. Also, as turbulence is a random process, there will be occasions when higher order aberrations are strong enough to further disrupt the laser beam profile.

Traditional Adaptive Optics (AO) Atmospheric Compensation

One technique that has been used to mitigate the fading problem in free space optical (FSO) laser communication links, is traditional adaptive optics (AO) atmospheric compensation. In principle, an AO compensation system on the transmitter can improve the condition of the beam at the receiver and help to reduce the fluctuation in received power. An adaptive optics compensation system on the receiver should also reduce the size of the spot at the detector plane.

Unfortunately, under conditions like those of the Wave-Optic Propagation Simulations described above, the performance of phase only AO compensation systems is very limited. The reason lies in the difference in the strength and distribution of the turbulence along a horizontal path compared to the vertical path of typical atmospheric AO compensation systems (e.g. astronomical adaptive optics). With turbulence roughly uniformly distributed along the path, phase errors close to the transmitter evolve into strong intensity variations at the far end of the path. Similarly, the light from the beacon on the receiver that is used as the source for the transmitter AO compensation system arrives at the transmitter with strong intensity variations. Portions of the received beam with very low intensity lead to erroneous phase measurements and corrupt the calculated wavefront. Even if the wavefront were measured perfectly and the correct phase compensation was applied to the outgoing laser beam, there would still be intensity variations across the laser beam at the receiver. This is because the transmitted laser beam is typically of uniform intensity while the laser beam that would produce uniform intensity at the receiver has intensity variations that mimic those of the laser beam from the beacon. This can be understood quite simply. Consider a portion of the laser beam from the beacon that is of low intensity. The light in that part of the laser beam has been diverged by the atmosphere to other parts of the beam. When the transmitted uniform beam passes through that same path, the light is converged by the atmosphere to produce a region of high intensity. What is needed to properly compensate the transmitter is full conjugation of the complex electric field of the laser beam from the beacon. While techniques exist for doing this, none have been demonstrated in the field.

At the receiver end, the compensation problem is not as daunting. There are the same problems due to regions of low intensity in the receiver pupil mentioned with regard to the transmit side of an AO compensation system. However, if the phase is correctly measured and compensated, then the laser beam spot at the detector should approach a diffraction limited size. This compensation system cannot do anything to reduce the size of variations in the power coming into the receiver. If, due to those variations, no light enters the receiver's entrance pupil, then no amount of phase compensation will increase the received signal power.

On top of these basic physical limitations of AO compensation systems for horizontal paths, there are also large and costly engineering problems due to the high temporal frequency characteristics of atmospheric turbulence. Measurements over a 16 km horizontal path show significant power in the full aperture tilt at temporal frequencies above 1 kHz. Higher order aberrations probably include even higher temporal frequencies. This requires an AO compensation system with a correction bandwidth of a least 1 kHz. Such response characteristics are beyond the performance of any existing AO systems.

Also, in general, some kind of reference wavefront is required for any type of wavefront measurement system, regardless of whether used as part of an OA subsystem in an FSO laser communication system, or some other kind of instrument or system. In any of the wavefront sensor configurations described above, the concept of the “reference wavefront” may be defined. In a phase-shifting interferometer, that reference wavefront is created by the physical surface of the reference flat or sphere used to form the interferometer. For a Hartmann wavefront sensor as described for example, in Applicant's U.S. Pat. No. 6,631,991, incorporated herein by reference as if set forth fully herein, the reference takes the form of the stored locations of the Hartmann spots for each subaperture—when the selected reference wavefront is input. The WFS control loop drives the spots toward these reference spot positions within each subaperture. In an ideal system, these subaperture reference positions would all be zero. In any real system, the null positions differ from zero. Typically, these positions are defined by causing a “perfect” plane wavefront to enter the WFS and recording the observed positions of the spots. This effectively calibrates the WFS reference positions so that the output of each subaperture is very close to zero in each axis when a plane wave is presented to the sensor. For many applications of wavefront sensing, the generation of this reference wavefront and its introduction into the wavefront sensor present significant difficulties.

Perhaps the most common problem to be overcome in creating a reference wavefront is that of non-common path optical aberrations. An example of this is the case is where the wavefront sensor uses one range of wavelengths to make its measurement, while it is desired to measure or correct a beam of a different wavelength. Because of chromatic aberrations in the optics, the wavefront that the sensor measures will differ from the wavefront that is to be corrected. In this case, there is no simple way to introduce the proper reference wavefront that will produce the correct measurement. While prior art laser communications systems provide good examples of this problem, the same problems are nevertheless present in all AO systems. Suppose the WFS has been calibrated with a plane wave and the spot reference positions recorded. This defines the null wavefront to be one that is flat at the plane of the WFS input. The control system will work to drive the wavefront at that plane toward the flat condition. Consider now the collection optic and receive fiber assembly. If the fiber is precisely at the infinity focus of the collection optic, a plane wave input will be optimally coupled into the fiber. If the fiber is slightly out of the correct focal plane, or if it is translated so that the plane wave that focuses onto the fiber is at an angle to the plane wave used to calibrate the sensor, then the control system will drive the wavefront to be plane at the WFS, but not optimally couple into the fiber. Further, if the beamsplitter between the WFS and the collection optic introduces different aberrations on the transmitted and reflected beams, then the control loop again will do the wrong thing by flattening the wavefront at the WFS.

Various several techniques also have been developed for addressing such “beam-fiber coupling” problems, which arise when the receiver fiber is slightly out of the correct focal plane or if it is translated so that the control system drives the wavefront to be planar at the WFS, but not optimally coupled into the fiber. For example, one compensation technique involves careful mechanical adjustment of the angle between the WFS reference and the boresight of the fiber so as to remove this error. Another compensation technique involves calibrating and applying aberrations of the beamsplitter to the WFS reference so that the null of the control loop produces the correct wavefront at the collection optic, rather than at the WFS entrance. However, both of these compensation techniques have the same shortcoming: the magnitude of the misalignment between the tracking sensor and the receiver is not fixed; it can vary with temperature, gravity loading, vibration, etc. Over the course of time, the alignment will degrade and the performance of the receiver system will be compromised. Optical aberrations are also dependent on mechanical alignment and will also vary with time.

Also, in a typical Hartmann wavefront sensor, the Hartmann spot positions are determined by using an algorithm such as “centroiding”. To avoid confusion with adjacent spots, it is necessary to extract from the image, a subarray that contains just one spot. This requirement limits the maximum allowable wavefront gradient in each subaperture of the Hartmann sensor. Furthermore, the gradient cannot be so large that the spot moves out of its subarray. This often leads to a limitation in overall system performance.

In summary, optical communications technology has rapidly advanced and is at the point where communication systems via fiber at 40 Gbs are commercially available. There are applications that prevent the use of a fiber connection, such as those involving moving platforms.

However, the transition of the light from fiber to free space presents problems. To achieve very high data rates requires the use of single mode fiber to prevent mode dispersion from corrupting the data stream. Conversion of the light signals to electrical signals at very high rates places constraints on the physical size of the electro-optical sensor. The capacitance of a detector is proportional to its area and, because of its inherent RC time constant, larger detectors are necessarily slower. For these reasons, to couple free space propagating light into a fiber or onto a detector requires that the light be concentrated into a very small region.

In the case of a single mode fiber, the mode diameter is similar to the diffraction limited spot size defined by the numerical aperture (NA) of the fiber. For efficient coupling, the NA of the coupling lens must match the NA of the fiber. This means that the focal spot needs to be close to diffraction limited to properly couple the light from free space to the fiber.

Unfortunately, propagation of light through the atmosphere introduces aberrations in the optical wavefront that prevent reaching the diffraction limit. To some extent, these aberrations can be corrected with adaptive optics. However, under conditions of strong turbulence (such as encountered in a free space link close to the ground), the disturbance of the optical beam include both phase and intensity components. Traditional phase only adaptive optics cannot compensate for these intensity disturbances. The result is that FSO links under these conditions experience deep signal fade events that lead to unacceptable interruptions in the communications link.

Therefore, there is a great need in the art for an improved method and apparatus for mitigating or compensating for the effects of laser carrier signal fluctuations (i.e. fading) detected at the receiver of FSO laser communication systems, while avoiding the shortcomings and drawbacks of prior art methodologies and apparatus.

OBJECTS AND SUMMARY OF THE PRESENT INVENTION

Accordingly, it is a primary object of the present invention to provide an improved method of and apparatus for mitigating or compensating for the effects of laser carrier signal fluctuations (i.e. fading) observed at the receiver of FSO laser communication systems, while avoiding the shortcomings and drawbacks of prior art methodologies and apparatus.

Another object of the present invention is to provide a free-space optical laser communication system supporting optically-separated signal transmission and reception channels, and employing laser beam speckle tracking mechanisms and speckle-to-fiber/detector locking mechanisms along the signal reception channels thereof for automatically stabilizing variations in the detected intensity of received laser beam carrier signals caused by atmospheric turbulence along said signal channels.

Another object of the present invention is to provide a free-space optical laser communication system having optically-separated signal transmission and reception channels at all terminal points in the communication system, wherein laser beam speckle tracking (i.e. following) mechanisms are employed in transmission channels of system to achieve a first level of optical signal intensity stabilization at signal detector of signal reception channel; and wherein speckle-to-fiber/detector locking mechanisms are employed along the signal reception channels of system to achieve a second level of optical signal intensity stabilization at the signal detector in the signal reception channel.

Another object of the present invention is to provide a method of and apparatus for automatically stabilizing variations in the intensity of received laser beam carrier signals caused by atmospheric turbulence along the signal reception channels of a free-space optical laser communication system supporting optically-separated signal transmission and reception channels.

Another object of the present invention is to provide a free-space optical laser communication system supporting optically-combined signal transmission and reception channels, and employing laser beam speckle tracking mechanisms and speckle-to-fiber/detector locking mechanisms along the signal reception channels thereof for automatically stabilizing variations in the intensity of received laser beam carrier signals caused by atmospheric turbulence along said signal channels.

Another object of the present invention is to provide a free-space optical laser communication system having optically-combined signal transmission and reception channels at all terminals in the communication system, wherein laser beam speckle tracking mechanism is employed in the signal transmission channels of all terminals in the communication system to achieve a first level of optical signal intensity stabilization at signal detector in each signal receiving channel, and wherein speckle-to-fiber/detector locking mechanism are employed in the signal reception channels of all terminals in the communication system to achieve a second level of optical signal intensity stabilization at signal detector in each signal reception channel.

Another object of the present invention is to provide a method of and apparatus for automatically stabilizing variations in the intensity of received laser beam carrier signals caused by atmospheric turbulence along the signal transmission channels of a free-space optical laser communication system supporting optically-combined signal transmission and reception channels.

Another object of the present invention is to provide a free-space adaptive optical laser communication system supporting optically-separated signal transmission and reception channels and employing laser beam speckle tracking mechanisms and speckle-to-fiber/detector locking mechanisms along the signal reception channels thereof for automatically stabilizing variations in the intensity of received laser beam carrier signals caused by atmospheric turbulence along said signal channels.

Another object of the present invention is to provide a free-space adaptive optical laser communication system having optically-separated signal transmission and reception channels at all terminals in the communication system, wherein wavefront sensing (WFS) and wavefront correction (WFC) mechanisms are employed in the signal transmission and reception channels of system (i.e. traditional adaptive optics) to improve the condition of the laser beam at the signal detector in the signal reception channel (i.e. reduce the size of the spot at the detector plane); wherein laser beam speckle tracking mechanisms are employed in the transmission channels of system so as to achieve a first level of optical signal intensity stabilization at signal detector in the signal reception channel; and wherein speckle-to-fiber/detector locking mechanisms are in signal reception channels of system to achieve a second level of optical signal intensity stabilization at signal detector in each signal reception channel.

Another object of the present invention is to provide a method of and apparatus for automatically stabilizing variations in the intensity of received laser beam carrier signals caused by atmospheric turbulence along the signal reception channels of a free-space adaptive optical laser communication system supporting optically-separated signal transmission and reception channels.

Another object of the present invention is to provide a free-space adaptive optical laser communication system supporting optically-combined signal transmission and reception channels and employing laser beam speckle tracking mechanism and speckle-to-fiber/detector locking mechanisms along the signal reception channels thereof for automatically stabilizing variations in the intensity of received laser beam carrier signals caused by atmospheric turbulence along said signal channels.

Another object of the present invention is to provide a free-space adaptive optical laser communication system having optically-combined signal transmission and reception channels at all terminals in the communication system, wherein wavefront sensing and wavefront correction mechanisms are employed in signal transmission and reception channels of all terminals in the communication system (i.e. adaptive optics) to improve the condition of the laser beam at the receiver (i.e. reduce the size of the spot a the detector plane); wherein laser beam speckle tracking mechanisms are employed in signal transmission channels of all terminals in the communication system to achieve a first level of optical signal intensity stabilization at signal detector of each receiving channel; and wherein speckle-to-fiber/detector locking mechanisms are employed in signal receiving channels of all terminals in the communication system to achieve a second level of optical signal intensity stabilization at signal detector of each receiving channel.

Another object of the present invention is to provide a free-space adaptive optical laser communication system supporting optically-combined signal transmission and reception channels and employing speckle-to-receiver-aperture locking mechanism along the signal transmission channels of the system and laser beam speckle tracking mechanisms and speckle-to-fiber/detector locking mechanisms along the signal reception channels thereof for automatically stabilizing variations in the intensity of received laser beam carrier signals caused by atmospheric turbulence along said signal channels.

Another object of the present invention is to provide a free-space adaptive optical laser communication system having optically-combined signal transmission and reception channels at all terminals in the communication system; wherein wavefront sensing and wavefront correction mechanisms are signal transmission and reception channels of all terminals in the communication system (i.e. adaptive optics) to improve the condition of the laser beam at the receiver (i.e. reduce the size of the spot a the detector plane); wherein speckle-to-receiver-aperture tracking mechanisms are employed in the transmission channel of the communication system and laser beam speckle tracking mechanism in the reception channels thereof, so as to achieve a first level of optical signal intensity stabilization at signal detector of each receiving channel; wherein speckle-to-fiber/detector locking mechanisms are employed in signal receiving channels of all terminals in the communication system so as to achieve a second level of optical signal intensity stabilization at signal detector of each receiving channel.

Another object of the present invention is to provide a free-space optical laser communication system, wherein the receiver to sends the communication link power signal that it is using to perform the Speckle-to-Fiber/Detector Locking Method back to the transmitter via the optical or other communications link, and the transmitter uses this signal as the figure of merit in an control algorithm that drives a spatial phase or intensity modulator operating in a control loop, so that a bright speckle in the transmitted laser beam is pointed and locked onto the receiver aperture, while the receiver loop locks a bright speckle in the received laser beam in its focal plane, onto the receive fiber (or detector).

Another object of the present invention is to provide a free-space optical laser communication system, wherein a compensation subsystem is provided that produces the highest possible resolution at the focal plane or maximizes the power into a communications fiber.

Another object of the present invention is to provide such a free-space optical laser communication system, wherein the WFS compensation subsystem employed therein defines a figure of merit that measures the performance of the compensation system, and that subsytem uses such a figure of merit, collected in real-time, as input to an algorithm that actively modifies the stored reference wavefront of the sensor (i.e. represented as reference spot positions on a Hartmann sensor) in such a way to maximize that figure of merit.

Another object of the present invention is to provide an improved wavefront sensing control method, wherein when a Hartmann spot has moved some threshold distance away from the center of a subaperture region, then the WFS control algorithm of the present invention automatically redefines the subaperture region so as to recenter the Hartmann spot with the redefined subaperture region.

Another object of the present invention is to provide an improved wavefront sensing control method, which is capable of producing the highest possible resolution at the focal plane or to maximize the laser beam power into a communications fiber.

Another object of the present invention is to provide an improved wavefront sensing control (WFS) method, wherein the WFS control method defines and computes, in real-time, a figure of merit that measures performance of the system and is provided, as input to an algorithm that (i) actively modifies the stored reference wavefront of the wavefront sensor (i.e. reference positions) in such a way to maximize that figure of merit, and also (ii) actively modifies the subarray location (i.e. subaperture position) in response to the motion of the spots.

Another object of the present invention is to provide such an improved wavefront sensing control (WFS) method, wherein the compensation process is carried out while the sensor and the correcting element (e.g. deformable mirror) are operating in a closed control loop.

Another object of the present invention is to provide such an improved wavefront sensing control (WFS) method, wherein all the reference spot positions are slightly translates while monitoring the power entering the communications fiber.

Another object of the present invention is to provide such an improved wavefront sensing control (WFS) method, wherein translation of the reference spots is equivalent to introducing full-aperture tilt into the reference.

Another object of the present invention is to provide such an improved wavefront sensing control (WFS) method, wherein by translating the spots to maximize the power, the method is eliminating the physical misalignment between the communications fiber and the wavefront sensor.

Another object of the present invention is to provide such an improved wavefront sensing control (WFS) method, wherein the optimization algorithm can be run continuously as a closed-loop control system, albeit at a lower bandwidth than the aberration correction loop.

Another object of the present invention is to provide such an improved wavefront sensing control (WFS) method, wherein a fast control loop and first and second slow control loops are provided, and wherein the action of the first slow loop is to continuously monitor the signal strength and to use this strength as a metric of the success of updates to the reference offsets in the WFS loop, and wherein the action of the second slow loop is to continuously monitor the signal strength and to use this strength as a metric of the success of updates to the subaperture locations in the WFS loop.

Another object of the present invention is to provide such an improved wavefront sensing control (WFS) method, wherein the second control loop actively and continuously modifies the subarray location (i.e. subaperture position) in response to the motion of the spots.

Another object of the present invention is to provide such an improved wavefront sensing control (WFS) method, wherein the second control loop averages the spot locations over several frames of data so as to determine which spots are threatening to move out of their subarrays (subaperture regions), and these subarray locations are then updated to bring the spot closer to the center of the subarray (subaperture region).

Another object of the present invention is to provide such an improved wavefront sensing control (WFS) method, wherein when employed in open-loop correction systems (continuously while the WFS sensor is operating), the WFS method can may be used to correctly measure and correct for large full aperture tilts induced in moving platforms b, updating the subaperture (subarray) locations.

Another object of the present invention is to provide such an improved wavefront sensing control (WFS) method, wherein when employed in closed-loop correction systems, the WFS method of active updating of the reference spot locations can address acquisition and pointing over a large field while maintaining control loop performance.

These and other objects of the present invention will become more apparently understood hereinafter and in the Claims to Invention appended hereto.

BRIEF DESCRIPTION OF THE DRAWINGS

For a more complete understanding of how to practice the Objects of the Present Invention, the following Detailed Description of the Illustrative Embodiments can be read in conjunction with the accompanying Drawings, briefly described below, wherein:

FIG. 1A is a schematic representation of an aircraft in communication with a space communication satellite by way of a free-space optical (FSO) laser communication system link, wherein various communication platform parameters are defined for use in modeling the atmospheric channel and the performance of the FSO laser communication system, including the altitude of the aircraft, the velocity of the aircraft, and the altitude and orbital path of the satellite, wherein the aircraft altitude determines the starting point of the communication link path, the aircraft velocity sets the rate at which the laser beam translates through the atmosphere, and the satellite altitude and path determine the apparent slew rate of the transmitted laser beam;

FIG. 1B is a schematic representation of a typical phase map used to calculate aberrations in a horizontal laser communications link propagating through strong atmospheric turbulence;

FIG. 1C is a time sequence of laser beam intensity distributions detected at the receiver of a laser communications satellite, calculated using a simple atmosphere model of propagation of a laser beam to the satellite at 1000 km altitude from an aircraft at 35,000 feet, wherein the bar in each frame in the time sequence represents 10 m at the satellite;

FIG. 1D is a typical Refractive Index Structure constant C_(n) ² profile generated by a simple atmosphere model, which is used to estimate the integrated turbulence of free space along a vertical path from the aircraft at various altitudes;

FIG. 1E is a plot of versus aircraft altitude, based on the results of one realization of turbulence estimates;

FIG. 1F1 is a plot of the atmospheric coherence length (Fried's r_(o) parameter) versus aircraft altitude, calculated using a wavelength of 1550 nm;

FIG. 1F2 is a plot of isoplanatic angle versus aircraft altitude, calculated using a wavelength of 1550 nm;

FIG. 1G is a graphical representation of a phase screen generated for an aircraft altitude at an elevation of 35,000 feet, calculated using the integrated C_(n) ² value of the link path;

FIG. 1H is the diffraction limited profile of the laser beam spot at the receiver of a satellite (left image) compared against a set of two examples of the laser beam spot profile at the receiver of a satellite, wherein each image covers about an area 100 m by 100 m, and wherein the laser beam profile in the center image shows that little power is lost due to wavefront aberration, while the laser beam profile in the right image shows that significant loss occurs due to dispersion of the laser beam;

FIG. 1I is a time series of laser beam power received at the receiver, calculated for a platform velocity of 200 m/s over about 0.1 s of flight time, showing at time 20 msec that there is one deep fade event of about 10 db contained in this simulated data set;

FIG. 1J is a calculated Zernike Distribution for the rms variation in each of the lowest 15 Zernike terms representing a computed compensation for phase aberrations generated in the laser beam as a result of laser beam pointing errors along the laser communication link, wherein the bulk of the aberration is shown contained in the first three Zernike terms (i.e. two Zernike terms required for tilt compensation, and one Zernike term required for defocus compensation);

FIG. 1K is a simulation of a typical laser beam spot intensity distribution pattern shown varying randomly at the plane of the receiver aperture as the turbulence along a 10 km communication link path changes through an atmosphere having C_(n) ²=1.0×10⁻¹⁴ m^(−2.3), and wherein the transmitter and receiver apertures are both 20 cm in diameter and the wavelength is 1.5 microns, and the full speckle image is 1 m×1 m;

FIG. 1L is a simulated spot of the received laser beam at the focal plane of the receiver, having the characteristics of a speckle image, wherein the scale of the (individual bright) speckle image is 1.25 μrad/pixel or 640 μrad for the full width of the image, and the diffraction limited spot size is 6 pixels, peak-to-null;

FIG. 1M is a simulated spot of the received laser beam at the focal plane of the receiver, showing a very dark portion of this speckle image falling onto the detector fiber core of the receiver, wherein the size and location of the detector or fiber core of the receiver is indicated as a small circle, and wherein (i) the fluctuations in the received intensity and (ii) the signal variations due to the moving speckle pattern at the focal plane, are not correlated;

FIG. 1N is a plot of log of the laser beam intensity (power) observed at the receiver detector versus the number of reiterations taken by the Laser Beam Speckle-Tracking algorithm hereof to find a new peak intensity to follow (at the detector), graphically illustrating that the laser beam intensity on the detector is monitored and when the intensity level falls below a predetermined threshold, then the line of sight of the laser beam is “shocked” and moved to a new location approximately a typical speckle spacing from its current location, and that this process is repeated until the detector intensity rises to an acceptable level, and thereafter the dithering of the received laser beam on the detector (i.e. fiber) is restarted, wherein one such “shocking” event is observed in FIG. 1N at about iteration number 220;

FIG. 1O is a plot of the probability density function (PDF) for the laser beam intensity data shown in FIG. 1N, clearly showing that the Laser Beam Speckle Tracking (i.e. Following) Mechanism of the present invention can effectively eliminate the long tail of low power fades observed in the raw (i.e. uncompensated) intensity data detected at the receiver, with occasional, intentional fades induced during the re-acquisition of speckle components (i.e. energy carrying beam spot components) during the Laser Beam Speckle Tracking Method of the present invention;

FIG. 1P is a set of calculated blur spots formed by introducing phase steps of 0.2, 0.5 and 1.0 waves in the entrance pupil of the communication receiver, illustrating that when a phase step is added to a random phase screen in the entrance pupil, the result is a modulation of the speckle pattern of laser beam spot detected by the fiber or detector;

FIG. 1Q is a series of speckle patterns observed within a detected laser beam spot, that are sinusoidally modulated at each location by adding phase steps of 0.2, 0.4, 0.6, 0.8 and 1.0 waves to a phase screen in the entrance pupil of the receiver, illustrating that the intensity of a given location within the speckle pattern (of the laser beam spot) can be maximized (i.e. optimized) by adjusting the size of the phase step, in accordance with the Speckle-to-Fiber/Detector Locking Method of the present invention carried out at the receiver of each terminal in the communication system;

FIG. 1R1 shows a plot of simulated log laser beam intensity (power) observed at the receiver detector versus the number of reiterations taken by the Speckle-to-Fiber/Detector Locking Mechanism of the present invention (using a 2×2 spatial phase modulation panel introduced in the entrance pupil of the receiver);

FIG. 1R2 shows a plot of simulated Log laser beam intensity (power) observed at the receiver detector versus the number of reiterations taken by the Speckle-to-Fiber/Detector Locking Mechanism of the present invention (using a 4×4 spatial phase modulation panel introduced in the entrance pupil of the receiver);

FIG. 1S is a plot of the probability density function (PDF) for the laser beam (speckle pattern) intensity data for simulations carried out using a 4×4 spatial phase modulation panel introduced in the entrance pupil of the receiver, showing that even when using a 4×4 (low spatial resolution) spatial phase modulation panel, having only sixteen (16) degrees of freedom, the Speckle-to-Fiber/Detector Locking Mechanism is capable of reducing significantly the occurrence of deep fades by virtue of its ability to lock the brightest laser beam speckle (in the laser beam spot image) onto the fiber/detector of the receiver, as turbulent atmospheric conditions move pas the entrance pupil of thereof;

FIG. 2A is a schematic representation of an aircraft carrying a multi-point FSO laser communication platform in accordance with the principles of the present invention, and supporting numerous bi-directional FSO laser communication links with multi-point FSO laser communication platforms carried by other aircrafts, terrestrial vehicles and naval carriers, shown arranged in accordance with a military command, control and communications (C3) environment;

FIG. 2B is a schematic representation of a city having buildings, towers, and other civil structures, on which are supported a multi-point FSO laser communication platforms in accordance with the principles of the present invention, each supporting numerous bi-directional FSO laser communication links with others multi-point FSO laser communication platforms carried by civil structures arranged within a civilian environment;

FIG. 3A is a schematic representation of a free-space optical (FSO) laser communication system supporting optically-separated signal transmission and reception channels, and employing Laser Beam Pointing Mechanisms along the signal transmission channels thereof for automatically pointing the transmitted laser beam towards the receiver aperture of the communication system, and Laser Beam Speckle Tracking Mechanisms and Speckle-to-Fiber/Detector Mechanisms along the signal reception channels thereof for automatically stabilizing variations in the detected intensity of received laser beam carrier signals caused by atmospheric turbulence along said signal channels;

FIG. 3B is a schematic optical diagram of the free-space optical (FSO) laser communication system of FIG. 3A, showing the transmitter and receiver module of each terminal in the system being arranged along optically-separated signal transmission and reception channels, and Laser Beam Pointing Mechanisms being employed along the signal transmission channels thereof for automatically pointing the transmitted laser beam towards the receiver aperture of the communication system, Laser Beam Speckle Tracking Mechanisms being employed in both the signal reception channels of system for automatically tracking (i.e. following) maximum intensity laser beam speckle and moving away from low intensity (i.e. black) laser beam speckle so as to achieve a first level of optical signal intensity stabilization at signal detector of the receiver modules, and also Speckle-to-Fiber/Detector Mechanism in the signal reception channels of system to automatically lock a maximum intensity laser beam speckle onto a fiber/detector, and thereby achieve a second level of optical signal intensity stabilization at the signal detector in the receiver modules of the system;

FIG. 3C is a schematic optical diagram of the transmitter module in each terminal of the laser communication system of FIG. 3A;

FIG. 3D is a schematic optical diagram of the receiver module in each terminal of the laser communication system of FIG. 3A;

FIG. 3E is a schematic representation of a 4×4 spatial phase modulation panel of used in the receiver module of each terminal in the laser communication system of FIG. 3A;

FIG. 3F is a schematic representation of a 4 cell (i.e. quadrant-type) signal detector used in the laser beam spot tracking mechanism in the receiver module of each terminal in the laser communication system of FIG. 3A;

FIG. 4A is a flow chart showing the steps involved in a first control method (i.e. algorithm) of generating a control signal for supply to the fast steering mirror (FSM) used in the Laser Beam Tracking Mechanism in the transmitter module of each the laser communication system of the present invention, wherein the values of offsets (K_(H) and K_(V)) used in the fast tracking control loop of the control algorithm are set at the time of system initialization;

FIG. 4A′ is a flow chart showing the steps involved in a second control method (i.e. algorithm) of generating a control signal for supply to the fast steering mirror (FSM) used in the Laser Beam Tracking Mechanism in the transmitter module of each the laser communication system of the present invention, wherein a slower control loop is provided within the algorithm's fast tracking control loop, that continuously and automatically optimizes the values of offsets (K_(H) and K_(V)) used in the control algorithm;

FIG. 4B is a flow chart showing the steps involved in third method of generating a proportional control signal for supply to the FSM used in the Laser Beam Tracking Mechanism in the transmitter module of each the laser communication system of the present invention, wherein the values of offsets (K_(H) and K_(V)) used in the fast tracking control loop of the control algorithm are set at the time of system initialization;

FIG. 4B′ is a flow chart showing the steps involved in fourth method of generating a proportional control signal for supply to the FSM used in the Laser Beam Tracking Mechanism in the transmitter module of each the laser communication system of the present invention, wherein a slower control loop is provided within the algorithm's fast tracking control loop, that continuously and automatically optimizes the values of offsets (K_(H) and K_(V)) used in the control algorithm;

FIG. 4C is a flow chart showing the steps involved in fifth method of generating a proportional plus integral plus derivative control signal for supply to the FSM used in the Laser Beam Tracking Mechanism in the transmitter module of each the laser communication system of the present invention, wherein the values of offsets (K_(H) and K_(V)) used in the fast tracking control loop of the control algorithm are set at the time of system initialization;

FIG. 4C′ is a flow chart showing the steps involved in sixth method of generating a proportional plus integral plus derivative control signal for supply to the FSM used in the Laser Beam Tracking Mechanism in the transmitter module of each the laser communication system of the present invention, wherein a slower control loop is provided within the algorithm's fast tracking control loop, that continuously and automatically optimizes the values of offsets (K_(H) and K_(V)) used in the control algorithm;

FIG. 4D is a flow chart showing the steps involved in a first method of generating a control signal for supply to the fast steering mirror (FSM) used in the Laser Beam Speckle Tracking Mechanism in the receiver module of each the laser communication system of the present invention, wherein the values of offsets (K_(H) and K_(V)) used in the fast tracking control loop of the control algorithm are set at the time of system initialization;

FIG. 4D′ is a flow chart showing the steps involved in a second method of generating a control signal for supply to the fast steering mirror (FSM) used in the Laser Beam Speckle Tracking Mechanism in the receiver module of each the laser communication system of the present invention, wherein a slower control loop is provided within the algorithm's fast tracking control loop, that continuously and automatically optimizes the values of offsets (K_(H) and K_(V)) used in the control algorithm;

FIG. 4E is a flow chart showing the steps involved in third method of generating a proportional control signal for supply to the FSM used in the Laser Beam Speckle Tracking Mechanism in the receiver module of each the laser communication system of the present invention, wherein the values of offsets (K_(H) and K_(V)) used in the fast tracking control loop of the control algorithm are set at the time of system initialization;

FIG. 4E′ is a flow chart showing the steps involved in fourth method of generating a proportional control signal for supply to the FSM used in the Laser Beam Speckle Tracking Mechanism in the receiver module of each the laser communication system of the present invention, wherein a slower control loop is provided within the algorithm's fast tracking control loop, that continuously and automatically optimizes the values of offsets (K_(H) and K_(V)) used in the control algorithm;

FIG. 4F is a flow chart showing the steps involved in fifth method of generating a proportional plus integral plus derivative control signal for supply to the FSM used in the Laser Beam Speckle Tracking Mechanism in the receiver module of each the laser communication system of the present invention wherein the values of offsets (K_(H) and K_(V)) used in the fast tracking control loop of the control algorithm are set at the time of system initialization;

FIG. 4F′ is a flow chart showing the steps involved in sixth method of generating a proportional plus integral plus derivative control signal for supply to the FSM used in the Laser Beam Speckle Tracking Mechanism in the receiver module of each the laser communication system of the present invention, wherein within the algorithm's fast tracking control loop, a slower control loop is provided that continuously and automatically optimizes the values of offsets (K_(H) and K_(V)) used in the control algorithm;

FIG. 5A is a flow chart showing the steps involved in a first illustrative embodiment of the method of present invention, for generating spatial phase modulation panel drive signals to be provided to the spatial phase modulation panel (i.e. deformable mirror) employed in Speckle-to-Fiber/Detector Locking Mechanism in the receiver module of each laser communication system of the present invention;

FIG. 5B is a flow chart showing the steps involved in a second illustrative embodiment of the method of present invention, for generating spatial phase modulation (SPM) panel drive signals to be provided to the spatial phase modulation panel (i.e. deformable mirror) employed in the Speckle-to-Fiber/Detector Locking Mechanism in the receiver module of each laser communication system of the present invention;

FIG. 5C is a flow chart showing the steps involved in a third illustrative embodiment of the method of present invention, for generating spatial phase modulation (SPM) drive signals to the supplied to the spatial phase modulation panel (i.e. deformable mirror) employed in the Speckle-To-Fiber/Detector Locking Mechanism in the receiver module of each laser communication system of present invention;

FIG. 5D is a flow chart showing the steps involved in a first illustrative embodiment of the method of present invention, for generating speckle-to-fiber/detector locking drive signals to the spatial phase modulation panel (i.e. deformable mirror) employed in the receiver module of each laser communication system of the present invention;

FIG. 6A is a flow chart showing the steps involved in a first illustrative embodiment of the method of present invention, for generating spatial intensity modulation (SIM) drive signals to be provided to the spatial intensity modulation panel employed in the Speckle-To-Fiber/Detector Locking Mechanism in the receiver module of each laser communication system of the present invention;

FIG. 6B is a flow chart showing the steps involved in a second illustrative embodiment of the method of present invention, for generating spatial intensity modulation (SIM) drive signals in the Speckle-To-Fiber/Detector Locking Mechanism in the receiver module of each laser communication system of the present invention;

FIG. 6C is a flow chart showing the steps involved in a third illustrative embodiment of the method of present invention, for generating spatial intensity modulation drive signals to be provided to the spatial intensity modulation panel employed in the Speckle-To-Fiber/Detector Locking Mechanism in the receiver module of each laser communication system of the present invention;

FIG. 6D is a flow chart showing the steps involved in a first illustrative embodiment of the method of present invention, for generating spatial intensity modulation drive signals to be provided to the spatial intensity modulation panel employed in the Speckle-To-Fiber/Detector Locking Mechanism in the receiver module of each laser communication system of the present invention;

FIG. 7A is a schematic representation of a free-space optical (FSO) laser communication system supporting optically-combined signal transmission and reception channels, and employing Laser Beam Pointing Mechanisms being employed along the signal transmission channels thereof for automatically pointing the transmitted laser beam towards the receiver aperture of the communication system, and Laser Beam Speckle Tracking Mechanism and Speckle-to-Fiber/Detector Locking Mechanisms employed along the signal reception channels thereof for automatically locking a maximum intensity laser beam speckle onto a fiber/detector, and thereby stabilizing variations in the detected intensity of received laser beam carrier signals caused by atmospheric turbulence along said signal channels;

FIG. 7B is a schematic optical diagram of the free-space optical (FSO) laser communication system of FIG. 37, showing the transceiver module of each terminal in the system being arranged along optically-combined signal transmission and reception channels, and Laser Beam Pointing Mechanisms being employed along the signal transmission channels thereof for automatically pointing the transmitted laser beam towards the receiver aperture of the communication system, Laser Beam Speckle Tracking Mechanisms being employed in both the signal reception channels of system to achieve a first level of optical signal intensity stabilization at signal detector of the transceiver modules, and also a Speckle-to-Fiber/Detector Mechanisms in the signal reception channels of system to automatically lock a maximum intensity laser beam speckle onto a fiber/detector, and thereby achieve a second level of optical signal intensity stabilization at the signal detector in the transceiver modules of the system;

FIG. 7C is a schematic optical diagram of the transceiver module in each terminal of the laser communication system of FIG. 3A;

FIG. 8A is a schematic representation of a free-space adaptive optical (FS-OA) laser communication system supporting optically-separated signal transmission and reception channels and employing Laser Beam Pointing Mechanisms along the signal transmission channels thereof for automatically pointing the transmitted laser beam towards the receiver aperture of the communication system, and Laser Beam Speckle Tracking Mechanisms and Speckle-to-Fiber/Detector Locking Mechanisms along the signal reception channels thereof for automatically locking a maximum intensity laser beam speckle onto a fiber/detector, thereby stabilizing variations in the intensity of received laser beam carrier signals caused by atmospheric turbulence along said signal channels;

FIG. 8B is a schematic optical diagram of the FS-AO laser communication system of FIG. 8A, showing the transmitter and receiver module of each terminal in the system being arranged along optically-separated signal transmission and reception channels having AO-compensation mechanisms, and Laser Beam Pointing Mechanisms being employed along the signal transmission channels thereof for automatically pointing the transmitted laser beam towards the receiver aperture of the communication system, Laser Beam Speckle Tracking Mechanisms being employed in the signal reception channels of system to achieve a first level of optical signal intensity stabilization at signal detector of the receiver modules, and also Speckle-to-Fiber/Detector Locking Mechanism in the signal reception channels of system to automatically lock a maximum intensity laser beam speckle onto a fiber/detector, and thereby achieve a second level of optical signal intensity stabilization at the signal detector in the receiver modules of the system;

FIG. 8C is a schematic optical diagram of the transmitter module in each terminal of the laser communication system of FIG. 8A;

FIG. 8D is a schematic optical diagram of the receiver module in each terminal of the laser communication system of FIG. 8A;

FIG. 8E is a flow chart illustrating the data processing steps carried out simultaneously along the plurality of a subapertures of the Shack-Hartmann Wavefront Sensor (WFS) according to an illustrative embodiment of the present invention, wherein a slow (low-bandwidth) control loop is provided to continuously monitor the subaperture signal strength and use this measured signal strength as a metric of success (i.e. performance) in optimally updating the reference position offsets within the fast WFS control loop;

FIG. 8F1 is a schematic representation of a Hartmann spot shown near the center of a subaperture region;

FIG. 8F2 is a schematic representation of a Hartmann spot shown moving away from the center of a subaperture region;

FIG. 8F3 is a schematic representation illustrating that when a Hartmann spot has moved some threshold distance away from the center of a subaperture region, then the WFS control algorithm of the present invention automatically redefines the subaperture region so as to recenter the Hartmann spot with the redefined subaperture region;

FIG. 9A is a schematic representation of a free-space adaptive optical (FS-OA) laser communication system supporting optically-combined signal transmission and reception channels and employing Laser Beam Pointing Mechanisms along the signal transmission channels thereof for automatically pointing the transmitted laser beam towards the receiver aperture of the communication system, and Laser Beam Speckle Tracking Mechanisms and Speckle-to-Fiber/Detector Mechanisms along the signal reception channels thereof for automatically locking a maximum intensity laser beam speckle onto a fiber/detector of the receiver, thereby stabilizing variations in the intensity of received laser beam carrier signals caused by atmospheric turbulence along said signal channels;

FIG. 9B is a schematic optical diagram of the FS-AO laser communication system of FIG. 9A, showing the transceiver module of each terminal in the system being arranged along optically-combined signal transmission and reception channels having AO-compensation mechanisms, and Laser Beam Pointing Mechanisms being employed along the signal transmission channels thereof for automatically pointing the transmitted laser beam towards the receiver aperture of the communication system, Laser Beam Speckle Tracking Mechanisms being employed in the signal reception channels of system to achieve a first level of optical signal intensity stabilization at signal detector of the transceiver modules, and also Speckle-to-Fiber/Detector Locking Mechanism in the signal reception channels of system to automatically lock a maximum intensity laser beam speckle onto the fiber/detector of the receiver, and achieve a second level of optical signal intensity stabilization at the signal detector in the transceiver modules of the system;

FIG. 9C is a schematic optical diagram of the transceiver module in each terminal of the laser communication system of FIG. 9A;

FIG. 9D is a schematic diagram of the FS-OA laser communication terminal used in the laser communication system of FIG. 9A;

FIG. 9E is a photograph of a tripod-mounted FS-OA laser communication terminal constructed in accordance with the architecture of FIG. 9B, wherein the transceiver telescope has a 20 cm aperture and is mounted beneath the optical bench of the terminal, and a high-spatial resolution AO (WFS/WFC) compensation subsystem is provided for communication link characterization/compensation, and a real-time subsystem is provided for stabilizing the intensity of laser beam carrier signals detected at the transceiver of the system;

FIG. 9F is a photograph of a compact FS-OA laser communication terminal also constructed in accordance with the architecture of FIG. 9B, wherein the transceiver telescope has a 15 cm aperture and is mounted beneath the optical bench of the terminal, and a high-spatial resolution AO (WFS/WFC) compensation subsystem is provided for communication link characterization/compensation, and a real-time subsystem having a Laser Beam Speckle Tracking Mechanism (including a fast steering mirror FSM)) and a Speckle-to-Fiber/Detector Locking Mechanism (including a deformable mirror as a spatial phase modulation panel) are provided for automatically locking a maximum intensity laser beam speckle onto a fiber/detector of the receive, thereby stabilizing the intensity of laser beam carrier signals detected at the transceiver of the system, in response to atmospheric turbulence;

FIG. 10A is a schematic representation of a free-space adaptive optical (FS-OA) laser communication system supporting optically-combined signal transmission and reception channels and employing Laser Beam Pointing Mechanisms along the signal transmission channels thereof for automatically pointing the centroid of the transmitted laser beam towards the receiver aperture of the communication system, Speckle-To-Receiver-Aperture Locking Mechanisms along the transmission channels for locking laser beam speckle in the transmitted laser beam onto the receiver aperture of the receiver module in the system, and Laser Beam Speckle Tracking Mechanisms and Speckle-to-Fiber/Detector Locking Mechanisms along the signal reception channels thereof for locking a maximum intensity laser beam speckle onto the fiber or detector of the receiver, thereby automatically stabilizing variations in the intensity of received laser beam carrier signals caused by atmospheric turbulence along said signal channels; and

FIG. 10B is a schematic optical diagram of the FS-AO laser communication system of FIG. 10A, showing the transceiver module of each terminal in the system being arranged along optically-separated signal transmission and reception channels having AO-compensation mechanisms, and Laser Beam Pointing Mechanisms being employed along the signal transmission channels thereof for automatically pointing the centroid of the transmitted laser beam towards the receiver aperture of the communication system, Speckle-to-Receiver-Aperture Locking Mechanisms being employed in the signal transmission channels of the system, and Laser Beam Speckle Tracking Mechanisms being employed in the signal reception channels of system to achieve a first level of optical signal intensity stabilization at signal detector of the transceiver modules, and also Speckle-To-Fiber/Detector Locking Mechanisms in the signal reception channels of system to lock a maximum intensity laser beam speckle onto the fiber or detector of the receiver, and thereby achieve a second level of optical signal intensity stabilization at the signal detector in the transceiver modules of the system.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS OF THE PRESENTS INVENTION

Referring to the figures in the accompanying Drawings, the various illustrative embodiments of the methods, systems and devices (i.e. instruments) of the present invention will be described in great detail, wherein like elements will be indicated using like reference numerals.

General Description of Laser Communications Link Fade Prevention Techniques of the Present Invention

Through research, Applicants have discovered that the spatial intensity variations in the received laser beam signal at the input plane of a receiver lead to temporal fluctuations in the total power received. FIG. 1K shows a simulation of a typical intensity distribution at the plane of the receiver aperture. This case is for a 10 km path through an atmosphere with C_(n) ²=1.0×10⁻¹⁴ m^(−2.3). The transmitter and receiver apertures are both 20 cm in diameter and the wavelength is 1.5 microns. The full image shown is 1 m×1 m. This intensity pattern varies randomly as the turbulence along the path changes.

Clearly, unless the receiver is made very large, there will be significant variations in the received power. The receiver acts to focus the received laser beam spot onto a detector, either directly or via a fiber optic. In either case, the focal spot (of the laser beam) must be very small to match the size of the detector or fiber core. Because of the phase aberrations in the input signal to the receiver, the focal spot is aberrated and will always be larger than the diffraction limit.

FIG. 1L shows the simulated spot of the received laser beam at the focal plane of the receiver. This simulated spot has the characteristics of a speckle image. The scale here is 1.25 μrad/pixel or 640 μrad for the full width of the image. The diffraction limited spot size is 6 pixels, peak-to-null. For details on the nature of speckle images produced by coherent laser beams, reference is made to WIPO Publication No. WO 02/43195 A2 by Applicants, which is incorporated herein by reference in its entirety.

Occasionally, a very dark portion of this speckle image will fall onto the fiber or detector of the receiver in the communication terminal. Such a case is shown in FIG. 1M. The small circle in FIG. 1M indicates the size and location of the detector or fiber core of the receiver. The fluctuations in the received intensity and the signal variations due to the moving speckle pattern at the focal plane are not correlated. Thus it will sometimes happen that at a time when the total received power is low a dark part of the speckle pattern will be over the detector. During such events the power reaching the detector may be reduced by a factor of more than 10³ relative to the average power. It is very difficult to design a communications link that can maintain a low bit error rate (BER) during these deep fade events.

A primary object of the present invention is to provide several techniques for the reducing signal fading in optical free space communication links. Fading of the signal due to the effects of atmospheric turbulence is a well-documented problem for the optical communications links, and has been studied in great detail at Adaptive Optics Associates, Inc., and described hereinabove.

In general, atmospheric turbulence induces wavefront aberrations that evolve so that both the phase and intensity distributions across the laser beam become non-uniform. This can lead to two problems for the communications link, namely: (1) the spatial intensity variations of the laser beam generated at the input plane of the receiver by channel turbulence lead to temporal fluctuations in the total power received at the detector; and (2) this spatial intensity pattern of the received laser beam varies randomly as turbulence along the laser beam path changes, producing a moving speckle image at the focal plane of the receiver, having a very dark portion which may fall upon the receiver's fiber or detector, causing the total received laser beam power to drop to extremely low levels when the dark part of the speckle pattern is passing over the detector. These problems will be described in greater detail below.

The new methods disclosed herein, while somewhat related to traditional AO systems, are based on a reassessment of the basic requirements of the free-space optical (FSO) laser communications system. AO based compensation systems are designed to produce a near diffraction limited spot after passage through a turbulent atmosphere. That design is driven by the requirement either for high resolution imaging (astronomical AO), or for maximum power density in the spot (directed energy weapons). If that could be achieved for a long horizontal path, it would satisfy the requirements for a stable communications link. However, communications link is neither an imaging system nor a weapon system. The requirement is not to minimize the spot size or maximize the power density but rather to minimize the variations in the power. In a typical communications link the average received power is more than adequate to achieve a very low BER. The power level is driven in part by the fading problem. To assure that the link will function during a fade to 5% of the average power, the average power level is set 20× the minimum needed for the required BER.

The methods described here are aimed at stabilizing the power on the detector rather than maximizing that power. The first method uses adjustment of the line of sight to try to keep a bright portion of the speckle pattern on the fiber or detector. The second method uses phase modulation in the pupil plane to modify the speckle pattern in ways that help to stabilize the received power.

Laser Beam Speckle Tracking (i.e. Following) Method and Mechanism of the Present Invention

The characteristics of the speckle pattern at the detector plane are determined by the atmospheric aberrations and the optical system that produced the pattern. The overall width of the pattern is related to the strength of the aberrations. The typical size of the individual speckles is similar to the size of the diffraction limited spot for particular optical system. Suppose that a speckle pattern such as the one shown in FIG. 1M falls on a fiber or detector. Assume that the pattern is significantly larger than the fiber or detector. If the laser beam speckle pattern is moved to three or more positions separated by distances less than a typical speckle size and the received power recorded at each location, then it is possible to calculate the local gradient of the intensity pattern. The laser beam speckle pattern is then moved in the direction of the positive gradient by a distance less than a speckle size and the estimation of the gradient repeated. If the speckle pattern did not change, then the result of repeating this process several times would be to move the speckle pattern so that a peak of a laser beam speckle was centered on the fiber or detector. In the realistic case of an evolving speckle pattern, as long as the gradient estimation and pattern repositioning is performed several times within the timescale of speckle evolution, the result will be to keep the fiber or detector close to the peak of the laser beam speckle. This objective can be achieved using essentially a hill climbing algorithm that uses dithering of the position of the laser beam peckle pattern to determine the direction towards the local maximum of intensity. Simulations of this algorithm were performed to assess its performance, as follows.

A single phase screen with a Kolmogorov spectrum was used in the simulation. A large screen was generated and subsections extracted along a linear path to simulate turbulence moving past the input aperture. Both the strength of the phase screen and the rate at which the screen is moved are input parameters to the simulation. For each extracted phase screen, the point spread function (PSF) was calculated. Dithering of the image was simulated by looking at the intensity at four nearby locations in the PSF. A plane was then fitted to the intensities at those points. This fit was then used to calculate the local gradient of the intensity. The effective line of sight of the system was then updated based on the gradient and a new phase screen extracted. This process was repeated typically for 1000 iterations. For each iteration, the power falling on the “fiber” or “detector” was recorded. At the same time, the power falling on an uncompensated detector using the same phase screens was calculated for comparison with the compensated power. FIG. 9G shows a typical result plotting the log of the power versus iteration number for the compensated (green) and uncompensated (red) data. This simulation used a 10 cm receiver at 1.5μ wavelength and a phase screen with Fried parameter (r₀) of 1 cm. The phase screen was moved by 1 mm between iterations, which is roughly equivalent, at a wind speed of 10 mph, to a time step of 225 μs.

These simulations uncovered a shortcoming of the Laser Beam Speckle Tracking (i.e. Following) Method of the present invention. The hill climbing algorithm was effective in keeping the peak of a laser beam speckle centered on the fiber or detector for considerable periods of time. However, it was determined that eventually a situation arises in which the laser beam speckle simply fades away without moving. The hill climbing algorithm is then cast adrift and must work its way up to the peak if a new laser beam speckle is encountered. If, in the course of following a speckle, the fiber or detector location had moved far from the average center of the PSF, then it may take a large number of iterations for the algorithm to find a new peak to follow.

A crude solution to this problem was implemented in the simulation, and used to generate the log intensity versus iteration number plot of FIG. 1N. The intensity on the detector was monitored and if it fell below some threshold, then the line of sight of the system was “shocked” and moved to a new location approximately a typical speckle spacing from its current location. This was repeated until the intensity rose to an acceptable level and then the dithering restarted. One such event is seen in FIG. 9G at about iteration 220. One further enhancement was added to the Method. If the distance of the detector from the average center of the PSF was above some value, then the fiber or detector was re-centered and the algorithm restarted. It should be noted that these “shock” events can lead to fades in the signal until a laser beam speckle is re-acquired. Fortunately, they are anticipated events and the communications link could be interrupted in an orderly fashion.

FIG. 1O shows the probability density function (PDF) for the data in FIG. 1N. This clearly shows that the Laser Beam Speckle Tracking Method can effectively eliminate the long tail of low power fades seen in the raw data. The price paid for this improvement are the occasional, intentional fades during speckle reacquisition.

The Laser Beam Speckle Tracking Method of the present invention is practiced in all the FSO laser communication systems shown in FIGS. 3A through 10B.

Speckle-to-Fiber/Detector Locking Method of the Present Invention

While the Laser Beam Speckle Tracking (i.e. Following) Method of the present invention showed clear improvement in the stability of received power, it was not entirely free of fading problems. An alternative approach was then examined using the same simulation tools, as described below.

It is well known that application of a phase step in the entrance pupil of an imaging system alters the diffraction limited image in a very particular way. FIG. 1P shows the PSFs resulting from the application of phase steps 0.2, 0.5 and 1.0 waves in the entrance pupil. If a phase step is added to a random phase screen in the pupil, the result is modulation of the speckle pattern. FIG. 1P shows the result of adding phase steps of 0.2, 0.4, 0.6, 0.8 and 1.0 waves to a phase screen that produces a speckle pattern. The result is that the intensity is modulated sinusoidally at each location. This phenomenon leads to the concept of the second approach to intensity stabilization. As shown in FIG. 1Q, the size of the phase step can be adjusted to maximize the intensity at a given location within a speckle pattern. As the speckle pattern evolves, the optimum phase step can be discovered using the same dithering and hill climbing technique used in the speckle following. Here, however, it is the size of the phase step that is dithered while monitoring the intensity on the detector. This very simple arrangement was simulated and was found to have the same problem that the speckle following technique exhibited. While it can maximize the brightness of a speckle on the detector, if that speckle fades away, then even its maximum brightness falls below acceptable levels.

The single variable phase step may be thought of as an adaptive optics (AO) system of sorts, having a single degree of freedom, but without the use of wavefront sensing. It is known that an AO system with a number of degrees of freedom larger than the number of coherence cells in the pupil can, in principle, compensate the image so that it is close to diffraction limit. A question which the present inventors asked is: Is there a system with a number of degrees of freedom between these limits that can stabilize the intensity of the laser beam sufficiently to be useful in a laser communications system? To answer that question, a simulation was constructed that used phase screens with 4, 9 and 16 degrees of freedom arranged as square arrays of 2×2, 3×3, and 4×4 independent piston levels in the pupil of the imaging system, respectively. The piston values of each element were dithered while the laser beam speckle intensity on the fiber or detector was monitored. Each element was then displaced in the direction that led to the highest power on the fiber or detector. The results of these simulations show that even with as few as four degrees of freedom it is possible to very significantly reduce the occurrence of deep fades. FIGS. 1R1 and 1R2 shows the log intensity plots for simulations of 2×2 and 4×4 phase modulators, respectively. The simulation parameters are the same as those for FIG. 1N. Thus, this Speckle-to-Fiber/Detector Locking Method described above appears to have great promise as an approach to the elimination of deep fades in a FSO laser communications link. Based on preliminary analysis of system performance, a 4×4 spatial phase modulator can keep the received power stabilized within a 10 db band under conditions that produce fades deeper that 30 db for an uncompensated system shown in FIG. 1S.

The Speckle-to-Fiber/Detector Locking Method described above is comparatively easy to implement using either spatial phase or intensity modulation at the entrance pupil of the receiver. The key element is the spatial phase or intensity modulator itself. There are several well known methods for introducing spatial phase or intensity modulation on a laser beam at the entrance pupil of the receiver. The most common are deformable mirrors (DMs) for spatial phase modulators, and liquid crystal modulators for spatial intensity modulators. Within the DM category, there are traditional piezoelectric driven mirrors and the more recently developed micro-machined (MEMS) versions. To determine which is best suited to this technique, a set of baseline requirements has been developed.

Employing Spatial Phase Modulation at the Entrance Pupil of the Receiver to Implement the Speckle-to-Fiber/Detector Locking Method of the Present Invention

Measurements of the temporal power spectrum of turbulence along horizontal paths often show significant power at frequencies above 1 kHz. To compensate for such turbulence, the correction system of the present invention (i.e. Speckle-to-Fiber/Detector Locking Mechanism) must have a response time at least ten times faster than the fluctuations in the atmosphere, or at most 100 μs. During one response time interval, it is necessary to perform all the dithering of the various degrees of freedom and the associated data processing.

For the case of the 4×4 spatial phase modulator this means that 16 cycles of the dither motion must occur in that interval. This implies a minimum frequency response of the spatial phase modulator of 160 kHz. Fortunately, the spatial phase modulator can be made quite small by positioning it at a relayed pupil image plane. Since the system is monochromatic, the phase modulator only needs a dynamic range of ±_ wave (the simulations described here limit the spatial phase modulator to this range).

There are currently available piezo-type actuators in sizes of several millimeters on a side that have resonant frequencies of 300 kHz and greater than 2 μm throw. By bonding thin, flat mirrors onto an array of these actuators, a suitable phase modulator can be fabricated. FIG. 3E shows a schematic view of a spatial phase modulator constructed using such components. The bonding of the mirrors may be done with the mirrors in contact with an optical flat to assure alignment.

Alternatively, there are developmental MEMS spatial phase modulators that approach 100 kHz resonant frequencies with a cell size of 0.3 mm. If such modulators prove practical, then they offer advantages in power consumption and cost.

Intensity data is required at rates similar to the dither frequency of the spatial phase modulator, but can easily be supplied by the communications link detector itself since it must be capable of operation at frequencies in the GHz range.

FIGS. 3D and 3E show an illustrative embodiment for implementing the Speckle-To-Fiber/Detector Locking Method of the present invention. As shown in FIG. 3D, it is optically similar to a traditional AO compensation system using a deformable mirror (DM) in a relayed pupil plane. It also includes a fast steering mirror (FSM) to remove the full aperture tilt component of the aberrations. This FSM control loop does not have to operate at the very high frequencies of the spatial phase modulator (DM), but it does need a control bandwidth approaching 1 kHz for operation under conditions of strong turbulence.

The data processing requirements for implementing the Speckle-to-Fiber/Detector Locking Method are modest. The intensity data steam does not need to be faster than about 1 MHz. During each dither cycle of the spatial phase modulator, the processor must test to see what direction of modulator motion produces a higher signal and apply an offset in that direction to the appropriate actuator. It also must assure that the offset signals to the phase modulator are kept within the ±½ wave region. This can be achieved by either a look-up table or a mathematical function that calculates the drive signal equivalent to one wave. Since modern processors can perform more than 1000 operations in the time for a single signal update, this is not seen as a problem.

Employing Spacial Phase Modulation at the Entrance Pupil of the Receiver to Implement the Speckle-to-Fiber/Detector Locking of the Present Invention

Another way of improving the information content of the laser beam speckle image received at the entrance pupil of an FSO laser communication receiver involves active masking of the entrance pupil of the receiver using spatial intensity modulation techniques.

A typical atmospherically induced wavefront has about 5 waves of aberration, peal-to-valley. During simulations, this wavefront can be generated in such a way that it has a Kolmogorov distribution. Examination of this wavefront reveals that there are regions within the wavefront that a relatively flat. These regions may, however, be tilted fairly significantly. If that region is extracted from the wavefront and its tilt removed, the residual wavefront error is found to be less than a ¼ wave. If the aperture of the imaging system were masked down to allow only this portion of the wavefront to enter the receiver's pupil aperture, then a much improved image would result.

It has been discovered that the masking of the pupil results in a marked improvement in laser beam image quality. However, it must be realized that the location and size of such masks depends on the instantaneous atmospheric wavefront in the pupil. For this technique to be practical, there must be a way to actively select the region to be allowed to form the image. One approach is to look for flat parts of the wavefront. Unfortunately, this would require knowledge of the shape of the wavefront, and there are difficulties associated with using wavefront sensing for this purpose. An alternative method is to use a mask that is dithered while the laser beam image quality is monitored at the entrance pupil. In accordance with the present invention, the entrance pupil can be broken up into many smaller regions and each of these alternately masked and unmasked. The laser beam spot image quality for each state of the mask is recorded and the mask for that region is set to the state that results in the higher image information content. This process is repeated for each of the sub-regions until the whole entrance pupil is covered. This process can be carried out using essentially a hill climbing algorithm with the goal of maximizing the laser beam spot image information.

To assess the performance of such a hill climbing algorithm, a computer simulation was constructed. The simulation models a 10 cm receiver aperture operating at a wavelength of 0.9 μm. Phase aberrations are introduced into the pupil of the receiver. These aberrations are generated by extracting parts of a large phase screen that has Kolmogorov distributed aberrations. Four separate regions are extracted and added to produce the pupil wavefront. The extraction points are moved across the large phase screen to simulate the effect of wind. Each of the four extraction regions moves in a different direction. Both the strength of the phase screen and the effective wind speed can be varied. The simulated pupil wavefront is then used to calculate the laser beam blur spot formed by the receiver telescope. This blur spot is then convolved with a high resolution image of the scene to simulate the effect of the aberration on laser beam spot image quality.

In this simulation, image contrast is used as a measure of laser beam image quality. A small (64×64 pixels) region of the laser beam image is extracted and the spatial gradient of the intensity is calculated. The metric of laser beam image quality is the ratio of the rms of the magnitude of the gradient divided by the average image intensity.

The intensity function in the receiver pupil is assumed uniform. The spatial light modulator (SLM) is simulated as an N×N array of square subapertures whose transmission can be set to 1.0 (open) or to some other value (closed). The closed value can be 0.0 to simulate a high contrast modulator or to larger values to simulate lower contrast modulators such as quantum-well devices.

The simulation begins with all subapertures in the open state. One at a time, each subaperture is closed and the change in image quality metric is recorded. If the metric drops, the subaperture is opened again, if the metric increases, then the subaperture is kept closed. After all the subapertures have been tested in this way, a new pupil wavefront is extracted from the large phase screen. The process is repeated. However for subsequent iterations, the state of each subperture is inverted from its current state and the effect on the metric is noted. Those changes that lead to increased metric are retained, all other subapertures are returned to their initial state.

Simulations have been run using modulators with resolutions from 5×5 to 20×20 with phase screens simulating turbulence with C_(n) ² in the range of 10⁻¹⁵ to 10⁻¹⁴ m^(−2/3). FIG. 3 shows the results of a simulation with a 7×7 element SLM and an atmospheric path with r_(o)=1.5 cm. The effective wind speed was such that the phase screens moved about 2.5 mm per iteration. For a wind speed of 10 mph this is equivalent to an iteration rate of 2 kHz. To simulate the laser beam sport image seen by a video rate sensor, the laser beam image data was averaged into bins of about 16 ms (60 Hz frame rate). The images are a time sequence of raw (bottom) and corrected (top) simulated video output. The image contrast metric for all 1000 iterations of this simulation is measured. It has been found that the image contrast metric is improved by an average factor of 2.44 with peak improvement of a factor of 4.6. If the image contrast is taken to be inversely proportional to the smallest image feature resolved, then an increase in contrast by a factor of 2.44 represents an increase in laser beam image information by a factor of 6.

Alternative pupil masking techniques and image quality metrics may be used. One approach that has been tested in the simulation is to use a pupil mask of fixed size whose position is dithered. Again, a hill climbing algorithm is used to determine the mask position that results in the best image quality. The performance of this technique appears to be similar to that of the subaperture mask. Several different image quality metrics have been tested in the simulation. These include the normalized rms amplitude of the gradient of the intensity, the maximum intensity gradient, and intensity gradient histogram based metrics.

The Speckle-to-Fiber/Detector Locking Method described above using spatial intensity modulation techniques at the entrance pupil of the receiver is comparatively easy to implement. The key element is the modulator itself. There are several well know methods for introducing intensity modulation on a beam. The most common are liquid crystal modulators and deformable mirrors (DMs). More recently, quantum well modulators have been developed that offer higher speed operation but with lower on-off contrast ratio. To determine which is best suited to this technique, a set of baseline requirements has been developed.

Measurements of the temporal power spectrum of turbulence along horizontal paths often show significant power out to frequencies well of 100 Hz. To compensate for such turbulence the correction system must have a response time at least ten times faster than the fluctuations in the atmosphere, or at most 1000 μs. During one response time interval it is necessary to perform all the dithering of the various degrees of freedom and the associated data processing. For the case of the 5×5 intensity modulator this means that 25 cycles of the dither motion must occur in that interval. This implies a minimum frequency response of the phase modulator of 25 kHz. Fortunately, the modulator can be made quite small by positioning it at a relayed pupil image plane.

There are currently available DM based spatial light modulators that are capable of switching at the rates required by this Speckle-to-Fiber/Detector Locking Technique. One example is the Texas Instruments DLP chip. This device uses micro-machining techniques to produce an array of small tilting mirrors. These mirrors may be flipped between the on state, in which they reflect the light into the remainder of the optical system, and the off state, in which the light is directed to an optical stop.

It has been found that at this time conventional liquid crystal modulators are too slow to be applied to this technique, being limited to speeds below 1 kHz. Quantum well modulators are capable of speeds well above 1 MHz, but provide contrast ratios of only 2:1 or 3:1 between states. However, initial simulation results indicate that the performance of this technique is only moderately degraded by using a 3:1 contrast ratio modulator.

The data processing requirements for using spatial intensity modulation techniques at the entrance pupil of the receiver to realize the Speckle-to-Fiber/Detector Locking Method of the present invention, are considerable, but within the capabilities of present PC processors. A 40×40 pixel high speed image sensor would have a data rate of 40 MHz. For each frame of data a image quality metric must be calculated. A 1 GHz processor can perform 25 operations on each pixel of the data. The algorithm used in the simulations performs about 5 operations per pixel. During each dither cycle of the modulator, the processor must test to see what state of the modulator produces a higher image metric and, on the basis of the signal change, set the state of the appropriate actuator. This part of the algorithm needs only to operate at the 25 kHz update rate and is not a problem.

These Speckle-to-Fiber/Detector Locking Methods of the present invention is practiced in all FSO laser communication systems shown in FIGS. 3A through 10B.

Speckle-to-Receiver-Aperture Locking Mithod of the Present Invention

Over long paths with significant turbulence, the laser beam spot pattern from a laser communication terminal may evolve into a speckle pattern. These spatial intensity variations at the input plane of the receiver lead to temporal fluctuations of the total power received. Applicant has constructed simulations of a typical intensity distribution at the plane of the receiver aperture (e.g. for a 10 km path through an atmosphere with C_(n) ²=1.0×10⁻¹⁴m^(−2/3)). In the simulation, the transmitter and receiver apertures are both 20 cm in diameter and the wavelength is 1.5 microns. The full image shown is 1 m×1 m.

The size of the individual speckles in the laser beam spot is about the diffraction limit of the full transmitter aperture, and the overall pattern size is set by the diffraction limit of the atmospheric coherence length. In general, an atmospheric compensation system with a single plane of correction (single DM) cannot correct both the phase and intensity in a beam that has passed through a long path with distributed turbulence. Therefore, even with compensation, it has been discovered that there will be intensity variations in the receiver input plane.

Under some circumstances, the laser beam speckle size will exceed the receiver aperture size. When this occurs, there is the possibility for significant signal fade events as a dark speckle covers the aperture. These deep fades adversely affect the communications link performance, because if no light enters the aperture, then no light will reach the fiber or detector at the receiver.

In situations where the laser beam has become a speckle pattern at the receiver, it is also true that the receiver is in the far-field of the transmitter aperture. This means that the transmitter appears as a point source to the receiver and that phase aberrations introduced close to the transmitter (including those of the transmitter AO system) have diffracted into a spherical wavefront with spatial intensity variations. Aberrations introduced onto beam close to the receiver will still be observed as phase variations across the laser beam, but since they have not evolved significantly they do not produce large intensity variations.

While the transmitter cannot affect the wavefront at the receiver, it can alter the speckle pattern at the receiver. That pattern is the result of the propagation of the wavefront from the transmitter aperture to the receiver. The phase of that wavefront is the addition of the atmospherically induced phase variations near the transmitter aperture and the wavefront generated by the transmitter. As has been demonstrated for the case of fade prevention at the receiver, it is possible to alter the speckle pattern considerably by changing the wavefront phase by less than a fraction of the wave. Thus another aspect of the present invention involves the communication of a figure of merit (typically the received power) from the receiver to the transmitter so that the transmitter can dynamically point and lock a maximum intensity speckle in. the transmitted laser beam onto the receiver aperture of receiver on the other side of the laser communication link, to optimize the transmission of laser beam power into the entrance pupil of the communication receiver. This method can be carried out by dithering the transmitted wavefront and monitoring a figure of merit used in the receiver to achieve fade prevention. Preferably, the figure of merit used by the algorithm is the received power at the other end of the laser communication link. This data may be sent to the transmitter via the optical communications link or by other communications means.

The situation is slightly more complicated if the transmitter is also operating as a receiver with its own fade prevention and/or AO system. If it is using just a fade prevention system, the algorithm will need to use both its local figure of merit and that sent from the receiver. Since there will typically be many degrees of freedom available in the spatial phase or intensity modulator, it will, in general, be possible to optimize the received power at both ends of the laser communication link.

If the transmitter is operating an adaptive optical correction system as part of its receiver, then the control algorithm will have to compromise between optimal local received power and fade elimination at the distant receiver. The phase modulation is required to effectively prevent fades at the far receiver as typically small enough that they will not reduce the local received power by a large factor. The algorithm can be tailored to achieve the lowest two-way BER for the full duplex link.

One issue to be considered is the latency associated with the transmission of the figure of merit data across the link. For a link length of 300 km the time from when the light leaves the transmitter and acquires its aberrations and when it reached the receiver is 1 ms. It will take another 1 ms for the data to return. To this must be added the latency of the figure of merit detection and calculation. This accumulated latency will limit the bandwidth of the fade prevention system.

This Speckle-to-Receiver-Aperture Locking Method of the present invention is practiced in the FSO laser communication system shown in FIGS. 10A and 10B.

Brief Summary of Speckle-Control Based Methods of Present Invention for Reducing Fading in FSO Laser Communication Systems

In summary, three generalized methods and apparatus have been described above for reducing fading of laser beam signal intensity in FSO laser communication systems, namely: (1) automatically tracking laser beam speckle at the receiver of any terminal in a FSO laser communication system so to avoid dark speckle spots from dwelling on the fiber or detector thereof and causing deep fades in the communication link; (2) automatically locking maximum intensity laser beam speckles onto the fiber or detector thereof so as to stabilize fluctuations in the intensity of a laser bean carrier detected at the receiver of the communication terminal of the present invention; and (3) automatically locking a maximum intensity laser beam speckle onto the receiver aperture of the laser communication system by measuring a figure of merit (e.g. signal intensity strength) at the receiver and communicating this information to the laser beam transmitter, so the transmitter can point a maximum intensity laser beam speckle onto the receiver aperture of the laser communication system and lock the same in a closed loop system. It is appropriate at this juncture to now describe five generalized embodiments of FSO laser communication systems shown in FIGS. in which such inventions can be embodied to provide remarkable levels of fade mitigation. Notably, each of these generalized system embodiments can be practiced in either military or civilian kinds of communication environments shown in FIGS. 2A and 2B, respectively.

Free-Space Optical (FSO) Laser Communication System Supporting Optically-Separated Signal Transmission and Reception Channels, and Employing Laser Beam Speckle Tracking Mechanisms and Speckle-to-Fiber/Detector Mechanisms Along the Signal Reception Channels Thereof for Automatically Stabilizing Variations in the Detected Intensity of Received Laser Beam Carrier Signals Caused by Atmospheric Turbulence Along Said Signal Channels

In FIG. 3A, there is shown a first illustrative embodiment of a free-space optical (FSO) laser communication system 1 in accordance with the principles of the present invention, wherein first and second communication terminals 2A and 2B are in communication by way of a super broad-band FSO laser beam communication link 3 having optically-separated signal transmission and reception channels, and wherein Laser Beam Speckle Tracking Mechanisms and Speckle-to-Fiber/Detector Mechanisms are employed along the signal reception channels of each communication terminal, for automatically stabilizing variations in the detected intensity of received laser beam carrier signals caused by atmospheric turbulence along said signal channels. As shown in FIG. 3A, each communication terminal has a transmitter (module) 4 and a receiver (module) 5, each of which will be described in detail below.

The transmitter in each terminal includes a telescopic transmitting aperture, and an optical train embedded with the components of a Laser Beam Pointing Mechanism realized by a fast steering mirror and a quad-cell detector of the kind shown in FIG. 3F, supporting optics, and a processor for carrying out one of the tracking algorithms shown in FIGS. 4A through 4C.

The receiver in each terminal includes a telescopic receiving aperture, and an optical train embedded with the components of two fade mitigation mechanisms of the present invention, namely: the Laser Beam Speckle Tracking Mechanism hereof realized using a FSM, a quad-cell detector shown in FIG. 3F, supporting optics, and a processor carrying out the tracking algorithms shown in FIGS. 4D through 4F; and the Speckle-To-Fiber/Detector Locking Mechanism hereof realized using spatial phase or intensity phase modulator shown in FIG. 3E, a receiving fiber, a single cell detector, and the processor carrying out (i) one of the spatial phase modulation (SPM) control signal generation algorithms shown in FIGS. 5A through 5D when using spatial phase modulation (SPM) techniques, or (ii) one of the spatial intensity modulation (SIM) control signal generation algorithms shown in FIGS. 6A through 6D when using spatial phase modulation (SPM) techniques.

The object of the Laser Beam Speckle Tracking Mechanism employed in both the signal reception channels of system is to automatically track or follow a maximum intensity laser beam speckle and move away from low intensity (i.e. black) laser beam speckles (that might fall onto the receiving fiber) so as to achieve a first level of optical signal intensity stabilization at signal detector of the receiver modules. The object of the Speckle-to-Fiber/Detector Mechanism in the signal reception channels of system is to lock a maximum intensity speckle in the received laser beam onto the receiving fiber, so as achieve a second level of optical signal intensity stabilization at the signal detector in the receiver modules of the system.

Notably, the purpose of the tracker/FSM control loop is different in the receiver and transmitter.

In the receiver, the tracker/FSM control loop is used to realized the Laser Beam Speckle Tracking Mechanism of the present invention which, when using the control algorithm of FIG. 4D, produces a signal that is proportional to the displacement of the received focal spot from the location of the receiver fiber or detector. This signal is typically generated by forming an image with a portion of the received light on a quad-cell or other position sensitive imaging device. The tracker (quad-cell detector) is calibrated so that its null position is optically identical to the position of the receiver fiber. The displacement signal is the input to a control loop algorithm that drives the steering mirror to keep the received focal spot on the fiber. Any of the standard control algorithms shown in FIGS. 4D through 4F may be used.

In the transmitter, the tracker/FSM control loop is used to realize the Laser Beam Pointing Mechanism. Here, the goal of the transmitter tracking control loop is to keep the transmitted beam on the receiver aperture. This is accomplished by having the tracker (quad cell detector) measure the displacement between the bore-sight direction of the transmitter and some optical signal received from the receiver. Any of the standard control algorithms shown in FIGS. 4A through 4C may be used. This optical signal can be a separate beacon, as shown in FIG. 3C, or it may be a portion of the light from the transmitter associated with the receiver. The hardware implementation of either tracker subsystem is the same, the only differences are the source of the light used by the tracker and the definition of the null position of the tracker.

Methods of Generating a Control Signal for Supply to the Fast Steering Mirror (FSM) Used in the Laser Beam Tracking/Pointing mechanism in the Transmitter Module of Each the Laser Communication System of the present Invention

In this simple form of a laser beam pointing/tracking system, a lens is used to form an image of the target (the beacon source or a portion of the communication link light) on a detector divided into four quadrants, as shown in FIG. 3F. Each quadrant produces a signal proportional to the intensity of light reaching the quadrant. The four resulting signals can be combined to produce two composite signals that are estimates of the position of the focal spot relative to the junction between the quadrants: Horizontal Position Estimator=(S1−S2)+(S4−S3)/(S1+S2+S3+S4)+K _(H) Vertical Position Estimator=(S1−S4)+(S2−S3)/(S1+S2+S3+S4)+K _(V)

Where K_(H) and K_(V) are positional offsets used to remove optical and mechanical alignment errors between the quad cell and the desired spot location.

If the light entering the tracking system has passed through a tilt control device, such as a steering mirror driven by piezoelectric actuators, these position estimators can be used as the error input signals to a control algorithm. Typical control system algorithms apply drive signals to the steering mirror that are proportional to a combination of the current error signal and the integral of previous error signals. The tilting action of the mirror provides optical feedback to the control system by altering the position of the image on the quadrant detector. Thus the system will tend to drive the focal spot to a location that equalizes the four signals, S1, S2, S3, S4, making the position estimator values close to zero. An offset (i.e. K_(H) and K_(V)) applied to the signals to cause the control loop to drive the spot to other locations. This allows the system to compensate for alignment differences between the tracker optical axis and the axis of the other portion of the optical system.

More complex arrangements are possible. A larger array of detector pixels may be used and their signals combined to form an estimate of the centroid of the intensity pattern of the focal spot. Alternatively, the signals may be used in a least-squares fitting algorithm or matched filter algorithm to provide estimates of the spot position.

If the focal spot has significant intensity structure (e.g. is a resolved image of the target) the signals from the detector may be used in a correlation algorithm. This algorithm performs a correlation of the instantaneous intensity distribution with an average or predefined target intensity distribution to form the estimate of target image location.

Definition and Maintenance of Tracking Sensor Boresight

In any of the tracker configurations shown above the concept of the “tracking sensor boresight” may be defined. In the case of the simple quadcell tracker, the position estimators are zero when the spot is centered on the junction between the four detector cells. A line that originates at this point and passes through the optical center of the imaging lens defines the null direction or boresight of the tracker system. Consider a simple optical system that includes a tracker and a fiber optic receiver such as might be used in a laser communications link. If the tracker and steering mirror are coupled in a closed loop control system, its action will be to drive the boresight of the tracker, as modified by reflection off the steering mirror, to point in the direction of the incoming light (shown by the red arrows in the figure). In general, the boresight of the tracker will not be in the same direction as boresight of the fiber receiver (defined by the line connecting the center of the fiber and the optical center of the collection optic. This misalignment of the two boresights can be corrected in several ways. Using optical metrology techniques, the two boresight directions may the mechanically adjusted into coincidence. Alternatively, the track sensing algorithm may include offsets (K_(H) and K_(V)) which can be modified to compensate for this misalignment. These offsets K are applied to the calculated position estimates and thus generally do not have a geometrical representation. In the “boresight tracking” case, the offsets are equivalent to an offset angle between the boresight of the tracker and the system boresight.

Each of these boresight tracking methods have the same shortcoming: the magnitude of the misalignment between the track sensor and the receiver is not fixed. It can vary with temperature, gravity loading, vibration, etc. Over the course of time the alignment will degrade and the performance of the receiver system will be compromised.

In the active updating tracking subsystem, this problem is addressed by enclosing the fast tracking loop in a slower loop that uses one of the hill-climbing algorithms previously discussed to continuously optimize the values of the offsets. In the next figure, a simple proportional control track loop is shown on the left and, on the right, a loop with the addition of the active update using a simple hill-climbing algorithm.

FIG. 4A describes the steps involved in a first method of generating a “proportional control” signal for supply to the fast steering mirror (FSM) used in the Laser Beam Tracking Mechanism in the transmitter module of each the laser communication system of the present invention. The control drive signal is simply proportional to the error signal. This is the simplest of the three algorithms. For stability reasons, the proportion or fraction of the error signal (k in the flow chart) must be less than one. For this reason, proportional control always has an offset error, that is, the control signal always is offset from the correct position. In this control algorithm, the values of offsets (K_(H) and K_(V)) used in the fast tracking control loop of the control algorithm are set at the time of system initialization.

In digital control systems such as those in the tracking loop, these are all implemented by applying a set of weights to the last N error measurements and summing the result with the current control drive signal (that is the integral part of the control). For example, if the weight of the current error measurement is k_(p), and all others are zero the algorithm is P plus I. To add derivative control, the weights for the last two error measurements would be (k_(p)+k_(d)) for the current one and −k_(d) for the error measurement form the last time step. This digital algorithm is very flexible and can implement many other forms of control and filtering of the error signal. Mathematically, it can be represented as D(t)=D(t ⁻¹)=G*(w ₀ *E(t ₀)+w ⁻¹ *E(t ⁻¹)+w ⁻² *E(t ⁻²)+. . . ) where D is the control loop drive signal, G is the overall loop gain, w_(n) is the weight associated with time step t_(n), and E(t_(n)) is the error signal from time step t_(n).

In FIG. 4A′, a second control method (i.e. algorithm) is described for generating a control signal for supply to the fast steering mirror (FSM) used in the Laser Beam Tracking Mechanism in the transmitter module of each the laser communication system of the present invention. In this control algorithm, a slower control loop is provided within the algorithm's fast tracking control loop, that continuously and automatically optimizes the values of offsets (K_(H) and K_(V)) used in the control algorithm.

FIG. 4B describes the steps involved in third method of generating a “proportional plus integral” control signal for supply to the FSM used in the Laser Beam Tracking Mechanism in the transmitter module of each the laser communication system of the present invention. Here, in addition to the proportional control signal, a term that is a fraction (k_(p)) of the integral of the past error signals is added. This has the effect of removing the offset error of proportional control alone. In this control algorithm, the values of offsets (K_(H) and K_(V)) used in the fast tracking control loop of the control algorithm are set at the time of system initialization.

In FIG. 4B′, a fourth method id described for generating a proportional control signal for supply to the FSM used in the Laser Beam Tracking Mechanism in the transmitter module of each the laser communication system of the present invention. In this control algorithm, a slower control loop is provided within the algorithm's fast tracking control loop, that continuously and automatically optimizes the values of offsets (K_(H) and K_(V)) used in the control algorithm.

FIG. 4C describes the steps involved in fifth method of generating a “proportional plus integral plus derivative” control signal for supply to the FSM used in the Laser Beam Tracking Mechanism in the transmitter module of each the laser communication system of the present invention. Here, another term may be added that is proportional to the time derivative of the error signal (multiplied by fraction k_(d) in the flow chart). The effect of this term is to provide anticipatory control, that is, the loop uses the current derivative to predict the error in the next time step. This can give higher loop bandwidth by reducing the latency, but increases the noise in the loop. In the transmitter, the goal is to keep the transmitted laser beam on the receiver aperture. So, the null or reference point in each of the above-described control algorithms, and the control loops which they maintain within the transmitter, is set such that the transmitted beam is collinear with the beacon laser beam sent from the receiver to enable laser beam pointing/tracking operations in the laser communication system. In this control algorithm, the values of offsets (K_(H) and K_(V)) used in the fast tracking control loop of the control algorithm are set at the time of system initialization.

In FIG. 4C′, a sixth method is described for generating a proportional plus integral plus derivative control signal for supply to the FSM used in the Laser Beam Tracking Mechanism in the transmitter module of each the laser communication system of the present invention. In this control algorithm, a slower control loop is provided within the algorithm's fast tracking control loop, that continuously and automatically optimizes the values of offsets (K_(H) and K_(V)) used in the control algorithm.

Methods of Generating a Control Signal for Supply to the Fast Steering Mirror (FSM) Used in the Laser Beam Speckle Tracking Mechanism in the Receiver Module of Each the Laser Communication System of the Present Invention

The control algorithms for the receiver shown in FIGS. 4D through 4F′, are similar to the basic control algorithms for the transmitter described in FIGS. 4A through 4C′, respectively, except for the “null point” (i.e. reference point) of the control system. In the transmitter, the goal is to keep the transmitted laser beam on the receiver aperture. So the null or reference point in each of the above-described control algorithms, and the control loops which they maintain within the transmitter, is set such that the transmitted beam is collinear with the beacon laser beam sent from the receiver to enable laser beam pointing/tracking operations in the laser communication system.

In FIG. 4D, a first method is described for generating a control signal for supply to the fast steering mirror (FSM) used in the Laser Beam Speckle Tracking Mechanism in the receiver module of each the laser communication system of the present invention. In this control algorithm, the values of offsets (K_(H) and K_(V)) used in the fast tracking control loop of the control algorithm are set at the time of system initialization.

In FIG. 4D′, a second method of generating a control signal for supply to the fast steering mirror (FSM) used in the Laser Beam Speckle Tracking Mechanism in the receiver module of each the laser communication system of the present invention. In this control algorithm, a slower control loop is provided within the algorithm's fast tracking control loop, that continuously and automatically optimizes the values of offsets (K_(H) and K_(V)) used in the control algorithm.

In FIG. 4E, a third method is described for generating a proportional control signal for supply to the FSM used in the Laser Beam Speckle Tracking Mechanism in the receiver module of each the laser communication system of the present invention. In this control algorithm, the values of offsets (K_(H) and K_(V)) used in the fast tracking control loop of the control algorithm are set at the time of system initialization.

In FIG. 4E′, a fourth method is described for generating a proportional control signal for supply to the FSM used in the Laser Beam Speckle Tracking Mechanism in the receiver module of each the laser communication system of the present invention. In this control algorithm, a slower control loop is provided within the algorithm's fast tracking control loop, that continuously and automatically optimizes the values of offsets (K_(H) and K_(V)) used in the control algorithm.

In FIG. 4F, a fifth method is described for generating a proportional plus integral plus derivative control signal for supply to the FSM used in the Laser Beam Speckle Tracking Mechanism in the receiver module of each the laser communication system of the present invention. In this control algorithm, the values of offsets (K_(H) and K_(V)) used in the fast tracking control loop of the control algorithm are set at the time of system initialization.

In FIG. 4F′, a sixth method is described for generating a proportional plus integral plus derivative control signal for supply to the FSM used in the Laser Beam Speckle Tracking Mechanism in the receiver module of each the laser communication system of the present invention. In this control algorithm, a slower control loop is provided within the algorithm's fast tracking control loop, a slower control loop is provided that continuously and automatically optimizes the values of offsets (K_(H) and K_(V)) used in the control algorithm.

Methods for Generating Spatial Phase Modulation Panel Drive Signals to be Provided to the Spatial Phase Modulation Panel (i.e. Deformable Mirror) Employed in Speckle-to-Fiber/Detector Locking Mechanism in the Receiver Module of each Laser Communication System of the Present Invention

The goal of the algorithms driving the Speckle-To-Fiber/Detector Locking Mechanism of the present invention is to find phase and/or intensity modulation pattern that, when added to the electric field at the entrance pupil of an optical system maximizes some metric function. For the case of optical communications, that metric is the time average communications signal level. In general, this is an optimization problem and many of the known optimization algorithms may be applied.

Define the Problem Variables:

-   E(x,y); The complex e-field at the entrance pupil -   M(x,y); Phase and/or intensity modulation applied to pupil -   S; Time averaged signal level

For the purposes of algorithm description it will first be assumed that E is constant in time. Issues related to real-time application of the algorithms will be addressed later. In the real system, S is “calculated” by allowing the collection optics to focus the electro-magnetic field onto the collection aperture. This may be simulated computationally by using standard point spread function (PSF) calculation methods to generate the spatial intensity pattern at the focal plane. This intensity distribution may then be integrated over the collection aperture. Typically, the PSF is calculated by taking the Fourier transform of the complex e-field at the aperture (E+M). The square of the amplitude of this transform then represents the intensity distribution.

In any realization of this system, M will be in the form an N×N array of subregions or subapertures within the full aperture. Each of these subregions represents a single element or pixel of the modulation device. Over each subregion the value of M is constant. For such an arrangement we then have N² independent degrees of freedom for the modulator and we must find the set of N² values (Mij) that maximize S.

First Illustrative Embodiment of the Method of the Present Invention, for Generating Spatial Phase Modulation Panel Drive to the Provided to the Spatial Modulation Panel (i.e. Deformable Mirror) Employed in Speckle-to-Fiber/Detector Locking Mechanism in the Receiver Module of Each Laser Communication System of the Present Invention.

FIG. 5A describes the steps involved in a first illustrative embodiment of the method of present invention, for generating spatial phase modulation panel drive signals to be provided to the spatial phase modulation panel (i.e. deformable mirror) employed in Speckle-to-Fiber/Detector Locking Mechanism in the receiver module of each laser communication system of the present invention. In this method, the hill climbing algorithm provides a simple approach to the above described optimization problem, wherein each degree of freedom is changed by a small amount in sequence and the resulting effect on S is measured. In this way the gradient of S with respect to the N² independent variables Mij is estimated. The values of Mij are then updated to move in the direction of the maximum upward gradient. This process is repeated until a set of Mij is arrived at, any change from which leads to a reduction in S. The algorithm has then found at least a local maximum of S.

Second Illustrative Embodiment of the Method of the Present Invention, for Generating Spatial Phase Modulation Panel Drive Signals to be Provided to the Spatial Phase Modulation Panel (i.e. Deformable Mirror) Employed in Speckle-to-Fiber/Detector Locking Mechanism in the Receiver Module of Each Laser Communication System of the Present Invention

FIG. 5B describes the steps involved in a second illustrative embodiment of the method of present invention, for generating spatial phase modulation panel drive signals to be provided to the spatial phase modulation panel (i.e. deformable mirror) employed in Speckle-to-Fiber/Detector Locking Mechanism in the receiver module of each laser communication system of the present invention. In this method, an alternative form of hill climbing algorithm may be used as well. In this case, rather than testing the response of S to each M_(ij) and then calculating their updated values, the value of each M_(ij) may be updated as soon as the response of S to that variable is determined.

Third Illustrative Embodiment of the Method of the Present Invention, for Generating Spatial Phase Modulation Panel Signals to the Supplied to the Spatial Modulation Panel (i.e. Deformable Mirror) Employed in Speckle-to-Fiber/Detector Locking Mechanism in the Receiver Module of Each Laser Communication System of Present Invention

FIG. 5C describes the steps involved in a third illustrative embodiment of the method of present invention, for generating spatial phase modulation drive signals to the supplied to the spatial phase modulation panel (i.e. deformable mirror) employed in the Speckle-To-Fiber/Detector Locking Mechanism in the receiver module of each laser communication system of present invention. In this method, a more sophisticated optimization algorithm, sometimes referred to as “stochastic parallel perturbative gradient ascent”, is used. The algorithm estimates the gradient of S with respect to M_(ij) in a somewhat different way. At each iteration, all of the values of M_(ij) are altered in a random way and the change in S is noted. If that change is represented by ΔS and the perturbation applied to M by ΔM_(ij), it follows that ΔS=S(M _(ij) +ΔM _(ij))−S(M _(ij)).   (1) To find the gradient of S with respect to each element of M (M_(mn)), a Taylor series expansion of Eq. 1 gives ΔS/ΔM _(mn) ≈∂S/∂M _(mn) +Σ∂S/∂M _(ij) *ΔM _(ij) /ΔM _(mn) , {i≠m, j≠n}.   (2)

If the perturbations are truly stochastic and random, the second term in eq. 2 averages to zero and thus ΔS/ΔM_(mn)≈∂S/∂M_(mn). Once again, once the gradient of S is estimated, the values of M_(ij) may be updated to move in the direction of maximum upward gradient.

Fourth Illustrative Embodiment of the Method of Present Invention, For Generating Speckle-To-Fiber/Detector Locking Drive Signals to the Spatial Phase Modulation Panel (i.e. Deformable Mirror) employed in the Receiver Module of Each Laser Communication System of the Present Invention

FIG. 5D describes the steps involved in a fourth illustrative embodiment of the method of present invention, for generating speckle-to-fiber/detector locking drive signals to the spatial phase modulation panel (i.e. deformable mirror) employed in the receiver module of each laser communication system of the present invention. This alternative method addresses a major limitation of simple gradient climbing algorithms, namely the potential for them to locate a local maximum of S, not the true global maximum. There are several ways to reduce the likelihood of the algorithm being trapped at a local maximum. One example is the simulated annealing algorithm. This algorithm is based on the process that occurs when a crystalline substance is slowly cooled through its freezing point, hence the name “simulated annealing”.

The flow of the “simulated annealing” algorithm is very similar to those described above. A random perturbation of the modulator values, M_(ij), is applied and the resulting change in S is noted. If S is increased, the perturbation is accepted unconditionally. However, if S is decreased, the perturbation may also be accepted, but with a probability given by P=Exp(ΔS/T),   (3) where T is a parameter of the algorithm. Because of the original relation of this algorithm to annealing, T is known as the temperature parameter. As the algorithm searches for the maximum and the value of S increases, the T parameter is slowly reduced. This lowers the probability of accepting downward changes in S. The effect of this finite probability of accepting negative changes in S is to provide occasional “shocks” to the system that act to throw the search point out of local maxima. At the beginning of the process, the search point may be far from the maximum and there may be many local maxima that could trap the algorithm. With a large T parameter, the probability of a shock is high, helping to prevent trapping. As the algorithm approaches the true maximum, the possibility of trapping is reduced, so the T parameter may be lowered.

The key to this process is the “annealing schedule”, the plan or algorithm used to determine when and by how much to reduce T. This depends intimately on the functional relation between S and M_(ij), which cannot easily be calculated since it also depends on the unknown, atmospherically affected E(x,y). Optimal annealing schedules must be developed via simulation or experiment.

First Illustrative Embodiment of the Method of the Present Invention, for Generating Spatial Intensity Modulation Drive Signals to be provided to the Spatial Intensity Modulation Panel Employed in Speckle-to-Fiber/Detector Locking Mechanism in the Receiver Module of Each Laser Communication System of the Present Invention

The basic Spatial Intensity Modulation (SIM) Control Signal Generation algorithms described in FIGS. 6A through 6D are similar to the SPM Control Signal Generation algorithms shown in FIGS. 5A through 5D and described above, respectively, except that these signals are calibrated for and supplied to spatially intensity modulation elements, rather than spatial phase modulation elements.

Free-Space Optical (FSO) Laser Communication System Supporting Optically-Combined Signal Transmission and Reception Channels, and Employing Laser Beam Speckle Tracking Mechanism and Speckle-to-Fiber/Detector Locking Mechanisms Along the Signal Reception Channels Thereof for Automatically Stabilizing Variations in the Detected Intensity of Received Laser Beam Carrier Signals Caused by Atmospheric Turbulence Along Said Signal Channels

In FIG. 7A, there is shown a second illustrative embodiment of a free-space optical (FSO) laser communication system in accordance with the principles of the present invention, wherein first and second communication terminals 11A and 11B are in communication by way of a super broad-band FSO laser beam communication link 12 having optically-combined signal transmission and reception channels, and wherein Laser Beam Speckle Tracking Mechanisms and Speckle-to-Fiber/Detector Mechanisms are employed along the signal reception channels of each communication terminal, for automatically stabilizing variations in the detected intensity of received laser beam carrier signals caused by atmospheric turbulence along said signal channels. As shown in FIG. 7A, each communication terminal has a transceiver (module) which will be described in detail below.

The transceiver in each terminal includes a telescopic transmitting/receiving aperture, and an optical train embedded with the components of a Laser Beam Pointing Mechanism realized by a fast steering mirror and a quad-cell detector of the kind shown in FIG. 3F, supporting optics, and a processor for carrying out one of the tracking algorithms shown in FIGS. 4A through 4C. The transceiver further includes an optical train embedded with the components of two fade mitigation mechanisms of the present invention, namely: the Laser Beam Speckle Tracking Mechanism hereof realized using a FSM, a quad-cell detector shown in FIG. 3F, supporting optics, and a processor carrying out the tracking algorithms shown in FIGS. 4D through 4F; and the Speckle-To-Fiber/Detector Locking Mechanism hereof realized using spatial phase or intensity phase modulator shown in FIG. 3E, a receiving fiber, a single cell detector, and the processor carrying out one of the spatial phase modulation (SPM) control signal generation algorithms shown in FIGS. 5A through 5D.

The object of the Laser Beam Speckle Tracking Mechanism employed in both the signal reception channels of system is to automatically track or follow a maximum intensity laser beam speckle and move away from low intensity (i.e. black) laser beam speckles (that might fall onto the receiving fiber) so as to achieve a first level of optical signal intensity stabilization at signal detector of the receiver modules. The object of the Speckle-to-Fiber/Detector Mechanism in the signal reception channels of system is to lock a maximum intensity speckle in the received laser beam onto the receiving fiber, so as achieve a second level of optical signal intensity stabilization at the signal detector in the receiver modules of the system.

In this particular illustrative embodiment of the FSO communication system of the present invention, the function of the tracker is now the same on each side. This is possible because the transmitter boresight and receiver boresight are the same (defined by the single fiber location). If the FSM keeps the received spot on the fiber, the transmitted beam will retrace the path of the received beam and reach the aperture on the other side.

Free-Space Adaptive Optical (FS-OA) Laser Communication System Supporting Optically-Separated Signal Transmission and Reception Channels and Employing Laser Beam Speckle Tracking Mechanisms and Speckle-to-Fiber/Detector Locking Mechanisms Along the Signal Reception Channels Thereof for Automatically Stabilizing Variations in the Intensity of Received Laser Beam Carrier Signals Caused by Atmospheric Turbulence Along Said Signal Channels

A laser communications terminal with an adaptive optic compensation system operating should act to maximize the power collected by the communications fiber. However, there may be occasions when the turbulence conditions are such that the compensation system performance degrades. This may occur when the correlation length in the atmosphere becomes less than the actuator spacing of the deformable device, meaning that the system does not have the spatial resolution to fully correct the wavefront, or when the strength of the aberrations exceed the dynamic range of the corrector. Whatever the cause, if the compensation system is operating with degraded performance it will no longer provide optimum communications signal power. The potential exists for very significant variations in the received signal (fades) because the focal spot will become aberrated.

In FIGS. 8A and 8B, there is shown a third illustrative embodiment of a free-space adaptive optical (FS-AO) laser communication system 20 in accordance with the principles of the present invention, wherein first and second communication terminals 21A and 21B are in communication by way of a super broad-band FSO laser beam communication link having optically-separated signal transmission and reception channels, and wherein a standard Adaptive Optics Subsystem and the Laser Beam Speckle Tracking Mechanism and Speckle-to-Fiber/Detector Mechanism of the present invention are employed along the signal reception channels of each communication terminal, for automatically stabilizing variations in the detected intensity of received laser beam carrier signals caused by atmospheric turbulence along said signal channels. As shown in FIG. 8A, each communication terminal 21A, 21B has a transmitter (module) 23 and a receiver (module) 24, each of which will be described in detail below.

The transmitter in each terminal includes a telescopic transmitting aperture, and an optical train embedded with the components of a standard Adaptive Optics Subsystem. and the Laser Beam Pointing Mechanism realized by a fast steering mirror and a quad-cell detector of the kind shown in FIG. 3F, supporting optics, and a processor for carrying out one of the tracking algorithms shown in FIGS. 4A through 4C′.

The receiver in each terminal includes a telescopic receiving aperture, and an optical train embedded with the components of a standard Adaptive Optics Subsystem, and the two fade mitigation mechanisms of the present invention, namely: the Laser Beam Speckle Tracking Mechanism hereof realized using a FSM, a quad-cell detector shown in FIG. 3F, supporting optics, and a processor carrying out the tracking algorithms shown in FIGS. 4D through 4F′; and the Speckle-To-Fiber/Detector Locking Mechanism hereof realized using spatial phase or intensity phase modulator shown in FIG. 3E, a receiving fiber, a single cell detector, and the processor carrying out (i) one of the spatial phase modulation (SPM) control signal generation algorithms shown in FIGS. 5A through 5D when using spatial phase modulation (SPM) techniques, or (ii) one of the spatial intensity modulation (SIM) control signal generation algorithms shown in FIGS. 6A through 6D when using spatial phase modulation (SPM) techniques.

The object of the Laser Beam Speckle Tracking Mechanism employed in both the signal reception channels of system is to automatically track or follow a maximum intensity laser beam speckle and move away from low intensity (i.e. black) laser beam speckles (that might fall onto the receiving fiber) so as to achieve a first level of optical signal intensity stabilization at signal detector of the receiver modules.

The object of the Speckle-to-Fiber/Detector Mechanism in the signal reception channels of system is to lock a maximum intensity speckle in the received laser beam onto the receiving fiber, so as achieve a second level of optical signal intensity stabilization at the signal detector in the receiver modules of the system.

The object of the Adaptive Optics (AO) subsystem employed in the FS-AO laser communication system of the present invention is to remove phase aberrations in the wavefront of the transmitted and received laser beam to that it is focused to its diffraction limit at the entrance pupil of the receiver of the communication terminal. In the illustrative embodiments, the AO subsystem employs a Shack-Hartmann WFS which may be thought of as an assemblage of trackers. Each subaperture generates a track signal for the its Hartmann spot. The measurement of the wavefront gradient is fully parallel. Subaperture-to-subaperture information transfer is not needed until the gradients are integrated or “reconstructed” to form an estimate of the optical path difference (OPD) or wavefront.

Overview of the Shack-Hartmann Wavefront Sensor (WFS) of the Present Invention Employing Active Updating of Reference Positions and Subaperture Positions on Said Wavefront Sensor

The WFS control method of the present invention is capable of producing the highest possible resolution at the focal plane or to maximize the laser beam power into a communications fiber. In each case, the WFS control technique defines and computes, in real-time, a figure of merit that measures performance of the system and is provided, as input to an algorithm that (i) actively modifies the stored reference wavefront of the wavefront sensor (i.e. reference positions) in such a way to maximize that figure of merit, and also (ii) actively modifies the subarray location (i.e. subaperture position) in response to the motion of the spots. Notably, this compensation process is carried out while the sensor and the correcting element (e.g. deformable mirror) are operating in a closed control loop.

A simple example of such a WFS control algorithm might be one that slightly translates all the reference spot positions while monitoring the power entering the communications fiber. By using, for example, a hill-climbing technique, the spot translation that maximizes the power may be found. Translation of the reference spots is equivalent to introducing full-aperture tilt into the reference. By translating the spots to maximize the power, the algorithm is eliminating the physical misalignment between the communications fiber and the wavefront sensor. Because this kind of misalignment may vary with temperature and mechanical loading on the optical system, this optimization algorithm be run continuously as a closed-loop control system, albeit at a lower bandwidth than the aberration correction loop.

Higher-order aberrations within the WFS system may also be eliminated by using a similar process. However, rather than translating the reference spots as a unit (tilt), each spot location would be modified by a different amount. For example, a focus shift in the reference would be produced by moving each spot toward or away from the center of the field by an amount proportional to its distance from that center. The key to this technique is using a system level “figure of merit” to actively and continuously modify the reference wavefront used by the sensor to achieve the best system performance. A variety of figures of merit may be used depending on the goal of the compensation system. Similarly, a variety of algorithms may be used to arrive at the reference spot locations that maximize the figure of merit.

The WFS control method of the present invention shown in FIG. 8E resolves both low and higher order aberrations by enclosing a fast WFS loop in a pair of (first and second) slower loops that uses one of the hill-climbing algorithms previously discussed to continuously optimize the values of the (wavefront) reference positions and subaperture locations on the Hartmann imager sensor. In the illustrative embodiment, the fast control loop is a simple proportional control loop, and the slow control loops use simple hill-climbing algorithms so as to actively update reference positions and subaperture locations, in real-time. Any combination of the tracking algorithms and hill-climbing algorithms may be used to implement this WSF control technique.

The action of the first slow loop the WFS control method is to continuously monitor the signal strength and to use this strength as a metric of the success of updates to the reference offsets in the WFS loop. If, for example, because of mechanical drift the boresight of the WFS gradually diverges from the boresight of the receiver, this loop should change the average WFS offsets (equivalent to full aperture tilt) to maintain alignment. If the position of the laser beam detection fiber drifts longitudinally, changing the focus, the reference positions should be altered by the slow loop in such a way that the wavefront they define has the correct defocus so that the light is still coupled into the fiber. It is important that the perturbations of the offsets used by the first slow loop be small so that they do not significantly degrade the performance of the fast tracking/control loop. Perturbations that produce 10-20% changes in the signal should be sufficient to yield a good measurement of the efficacy of a given perturbation.

The action of the second slow loop is to actively modify the subarray location (i.e. subaperture position) in response to the motion of the spots. The average spot locations over several frames of data can be used to determine which spots are threatening to move out of their subarrays (subapertures). These subarray (subaperture) locations can then be updated to bring the spot closer to the center of the subarray (subaperture). This process is used continuously while the sensor is operating. In a open-loop wavefront sensing application, this technique may be used to deal with large but slowly varying aberrations. An example might be wavefront sensing from a moving platform, as shown in FIG. 1A. The platform motion may induce large full aperture tilts that would exceed the dynamic range of the sensor. By updating the subarray locations, this large tilt may be correctly measured. In closed-loop correction systems, this technique, coupled with the active updating (i.e. modification) of the reference spot locations described above, may address acquisition and pointing over a large field while maintaining control loop performance.

Description of Plural Data Processing Operations Along the Subapertures of the Shack-Hartmann Wavefront Sensor (WFS) of the Present Invention Employing Active Updating of Reference Positions and Subaperture Positions on Said Wavefront Sensor

Referring to FIG. 8E, the data processing steps carried out simultaneously along the plurality of a subapertures of the Shack-Hartmann Wavefront Sensor (WFS) will now be described in accordance with the principles of the present invention. Details concerning the structure and function of conventional Shack-Hartmann (Hartmann) wavefront sensors are disclosed in detail in Applicant's prior U.S. Pat. No. 6,631,991, supra.

As shown in FIG. 8E, the WFS control process of the present invention for use with the WFS employed in the laser communication systems disclosed herein comprises two control loops, namely: a fast control loop, and a slow control loop. The fast loop performs basic WFS data processing operations across the Hartmann wavefront sensor, including reconstruction of the wavefront (OPD function) and generation of wavefront error correction signals delivered to the wavefront correction device (e.g. deformable mirror array). The slow (low-bandwidth) control loop continuously monitors the subaperture signal strength and uses this measured signal strength as a metric of success (i.e. performance) in optimally updating the reference position offsets and location, within the fast WFS control loop.

Description of Data Processing Operations Carried Out Within the Fast Control Loop of the WFS Control Method of the Present Invention

Initializing System and WFS Control Algorithm/Method

The first step before entering the Fast Control Loop of the WF Sensing Control Process shown in FIG. 8E, is to initialize the wavefront sensing control process (or subsystem). Then, the reference positions are set to initial values, and to subsequently modified (i.e. updated by the slow control loop of the WFS control process.

Capturing Image

Image data containing Hartmann spot images is then collected. Offset and gain correction, based on dark and flat field images, is then performed.

Estimating Hartmann Spot Positions

Estimating the position of spots is a basic measurement performed in the Shack-Hartmann sensor. Centroid calculation is the simplest technique to perform this measurement. According to the Centroid technique, for an array of detector pixels centered on a spot, weights proportional to the X and Y distance of the pixel from the center of the array are applied to the signal of that pixel. Then, the weighted sum of signals are divided by the unweighted sum of signals so as to produce a position estimate. This technique is analogous to a center of mass calculation. P _(x) =ΣS _(i,j) Wx _(i,j) /ΣS _(i,j) where P_(x) is the x position estimate,

-   S_(ij) is the signal from the i,j^(th) pixels, and -   Wx_(i,j) is the x weight of the i,j^(th) pixel, which is equal to     the x coordinate of the pixel relative to the center of the array;     Wx_(i,j)=i−i₀, where i₀ is the index of the central pixel.

More sophisticated techniques can be used to improve the accuracy of the spot position estimate. An example is “spot fitting”. A least squares fit of a suitable function to the signal values in the pixel array is performed. An example of such a function is F(x,y)=Ax ² +Bxy+Cy ² which approximates the shape of a typical spot around it's peak. Once the values of A, B, and C are found from the fitting operation, the location of the maximum of F is calculated and used as the spot position estimate. This approach is much less sensitive to errors in offset subtraction and detector noise than the technique of centroiding, described above. Position repeatability of better than 1/100 pixel can easily be achieved with the method of the present invention. Subtracting Reference Spot Positions (“Offsets”)

Hartmann wavefront sensing is a differential measurement. It is necessary to supply a set of reference spot positions that represent the null or “perfect” wavefront. This is typically done by injecting a flat wavefront directly into the sensor and recording the spot positions. These spot positions are known as the reference positions (i.e. “offsets”). To make a measurement of a wavefront, it is necessary to determine the displacement of the spots from their reference positions (or offsets), hence the subtraction of the reference positions. Spot motion is proportional to wavefront gradient.

In the case of active updating of reference positions, the values of these positions are altered while making measurements. The goal is to optimize the reference positions based on some metric of system performance. This addresses the possibility that the “perfect” wavefront is not the optimal reference because of opto-mechanical manufacturing errors and thermal and gravitational drift.

Reconstructing the Wavefront (OPD Function) by Integrating Gradient Values Among WFS Subapertures

The differential spot motions are proportional to the gradients of the wavefront (i.e. Spot motions ∝∇φ). To determine the shape of the wavefront, it is necessary to perform some form of numerical integration. The wavefront gradient equation ∇φ is expressed mathematically as follows:

${\nabla{\phi\left( {x,y,z} \right)}} = {{\frac{\partial\phi}{\partial x}\hat{x}} + {\frac{\partial\phi}{\partial y}\hat{y}} + {\frac{\partial\phi}{\partial z}\hat{z}}}$ There are many techniques for solving this equation for φ. The most common is the matrix reconstructor. A matrix is generated that relates the gradients to the wavefront. This matrix simply forms the differences between adjacent wavefront values to represent the gradients. G=M*φ, where G is the vector of gradients, M is the matrix, and φ are the wavefront values. If each side of the equation is multiplied by M⁻¹, the following expression is obtained: M ⁻¹ *G=M ⁻¹ *Mφ=φ. Thus, M⁻¹ is the matrix that, when multiplying G, produces the desired wavefront values, φ. This “reconstructor matrix” is pre-generated. Other techniques for solving the wavefront gradient equation ∇φ for φ include least-squares fitting, maximum entropy, and iterative relaxation approaches. Updating DM Drive Signals

Once the wavefront values are estimated, they are used to calculate the updates to the deformable mirror drive signals. In a closed loop system, the measured wavefront represents the residual error signal for the control loop. It is the difference between the wavefront and the waveform that is represented by the reference positions. This error signal is the input to any of the control algorithms (proportional, PI, PID, etc.)

Calculating the Performance Metric

After updating the DM drive signals, the system performance metric S ( e.g. based on lasercom signal level, Power on target, Image sharpness) is calculated.

Terminating the AO Loop

The WFS control process then determines whether or not to terminate the AO loop. If the control process is not to be terminated, then the control process within the fast loop returns to capturing the next Hartmann spot image, and proceeding forward along the data processing path of the fast control loop, as shown in FIG. 8E. If the control process is to be terminated, then the control process within the fast loop is terminated and the control process is finished, as shown in FIG. 8E.

Description of Data Processing Operations Carried Out Within the First Slow Control Loop of the WFS Control method of the Present Invention

Calculating Temporal Average of the System Performance Metric

After updating the DM drive signals, the system performance metric S is calculated, then the sloop control loop is entered, by calculating the temporal average of the performance metric S. Such temporal averaging can be carried out using a low-pass data filter.

Calculating the Gradient of the Performance Metric with Respect to Each Component of the Reference Position Gradient

The gradient of the performance metric S is calculated with respect to each component of the reference positions, namely: G _(i,j) =∂S/∂Pref _(i,j) Updating Reference Positions

The reference positions are then updated as follows: Pref′ _(i,j) =Pref _(i,j) +K G _(i,j) ∂Pref _(i,j) Perturbing the Reference Positions

Then, the reference positions are perturbed as follows: Pref ^(P) _(i,j) =Pref ^(|) _(i,j) +∂Pref ^(|) _(i,j) Thereafter, these perturbed reference position values are used to set the reference position values in the fast control loop. Description of Data Processing Operations Carried Out Within the Second Slow Control Loop of the WFS Control Method of the Present Invention

The second slow (low-bandwidth) control loop, shown on the right side of the flow chart of FIG. 8E, addresses problem of large drifts in spot positions that might move them out of their defined subaperture regions. The details of this second slow control can be understood with reference to FIGS. 8E through 8F3.

In general, a Shack Hartmann wavefront sensor (WFS) produces an array of spot images on an image detector. Subregions of the image detector are defined around each spot. These are defined by an array of 4-element vectors. Each element contains the pixel coordinates of the lower left corner of the subarray and its X and Y dimensions.

Thus, for the subarray shown in FIG. 8F1, the entry in the list of subapertures would be (0, 0, 5, 5) indicating that the subarray starts at coordinate 0,0 and is 5×5 pixels. Every time the location of this spot is to be measured, this subarray of pixels is extracted and the spot position algorithm used to estimate the spot location. Ideally, the spot location should be close to the center of the subarray of pixels. In this case shown in FIG. 8F1, the center of the subarray, S_(cent), is at (2.5, 2.5). The spot is shown at the center as well so the spot position S_(pos) would be (2.5, 2.5).

In the presence of wavefront distortions, the spot moves away from its null position. In general, the distance of the spot from its null position may be divided into two portions based on the temporal characteristics of the distortion causing the motion. There is a high temporal frequency component that is typically due to atmospherically induced distortion. In addition, there may be a slowly varying component due to thermal and mechanical drift in the optical system. Thus these portion of the spot can be represented as follows: S_(pos)(t)=S_(atm)(t)+S_(drift)(t)

If a time average of S_(pos) is formed over a time long compared to the frequency content of S_(atm), then the contribution of S_(atm) will tend to zero average. Time Avg.(S_(pos))≈S_(drift). In FIG. 8F2, S_(pos) might be (3.5, 3.5) while S_(cent) has not changed.

As shown, the first data processing step in the second slow loop stems from the fast control loop where spot positions are estimated. If it is found that the time average position of the spot has moved far from the center of the subarray, then the location of the subarray can be updated: If [(S_(drift)−S_(cent))>Threshold] then [S_(cent)=S_(drift)]. This is implemented by altering the list of subarray locations to best center the subarrays on their respective time average spot locations. The resulting image would then be that shown in FIG. 8F3, and S_(pos) is still (3.5, 3.5) but S_(cent) has been updated to (3.5, 3.5). For typical systems, the value of Threshold may be in the range of 3-5 pixels and the subarray list would be updated at a rate of about once per second. These values can be adjusted over a wide range to meet the particular requirements of a given system.

Then, the second slow control loop involves comparing the Time Average Spot Positions S_(pos) to the Subaperture Center Locations S_(cent), and then returning to the fast control loop where subaperture locations are set with currently optimal values in accordance with the principles of the present invention.

By continuously updating the positions of those subaperture regions in accordance with the present invention, the dynamic range of the spot position measurement is significantly increased. In addition, the performance of the spot location estimation algorithm is best when all the light of the spot is captured. With a spot near the edge of the subaperture region, some of the light will spill over into the adjacent subaperture leading to errors in both measurements. By virtue of this second control loop in the WFS control algorithm of the present invention, subaperture locations are updated slowly so that fast changes in the wavefront due to turbulence, etc. do not affect the locations but slow drifts in spot location lead to adjustment of the subapertures.

Free-Space Adaptive Optical (FS-OA) Laser Communication System Supporting Optically-Combined Signal Transmission and Reception Channels and Employing Laser Beam Speckle Tracking Mechanisms and Speckle-to-Fiber/Detector Mechanisms Along the Signal Reception Channels Thereof for Automatically Stabilizing Variations in the Intensity of Received Laser Beam Carrier Signals Caused by Atmospheric Turbulence Along Said Signal Channels

In FIGS. 9A and 9B, there is shown a fourth illustrative embodiment of a free-space adaptive optical (FS-AO) laser communication system 30 in accordance with the principles of the present invention, wherein first and second communication terminals 31A and 31B are in communication by way of a super broad-band FSO laser beam communication link 32 having optically-combined signal transmission and reception channels, and wherein a standard Adaptive Optics Subsystem and the Laser Beam Speckle Tracking Mechanism and Speckle-to-Fiber/Detector Mechanism of the present invention are employed along the signal reception channels of each communication terminal, for automatically stabilizing variations in the detected intensity of received laser beam carrier signals caused by atmospheric turbulence along said signal channels. As shown in FIGS. 9A and 9B, each communication terminal has a transceiver (module) 31A which will be described in detail below.

The transceiver in each terminal includes a telescopic transmitting/receiving aperture, and an optical train embedded with the components of a standard Adaptive Optics Subsystem. and the Laser Beam Pointing Mechanism realized by a fast steering mirror and a quad-cell detector of the kind shown in FIG. 3F, supporting optics, and a processor for carrying out one of the tracking algorithms shown in FIGS. 4A through 4C. The optical train of the transceiver further includes the components of the two fade mitigation mechanisms of the present invention, namely: the Laser Beam Speckle Tracking Mechanism hereof realized using a FSM, a quad-cell detector shown in FIG. 3F, supporting optics, and a processor carrying out the tracking algorithms shown in FIGS. 4D through 4F; and the Speckle-To-Fiber/Detector Locking Mechanism hereof realized using spatial phase or intensity phase modulator shown in FIG. 3E, a receiving fiber, a single cell detector, and the processor carrying out (i) one of the spatial phase modulation (SPM) control signal generation algorithms shown in FIGS. 5A through 5D when using spatial phase modulation (SPM) techniques, or (ii) one of the spatial intensity modulation (SIM) control signal generation algorithms shown in FIGS. 6A through 6D when using spatial phase modulation (SPM) techniques.

The object of the standard Adaptive Optics Subsystem is to remove phase aberrations in the wavefront of the transmitted and received laser beam to that it is focused to its diffraction limit at the entrance pupil of the receiver of the communication terminal.

The object of the Laser Beam Speckle Tracking Mechanism employed in both the signal reception channels of system is to automatically track or follow a maximum intensity laser beam speckle and move away from low intensity (i.e. black) laser beam speckles (that might fall onto the receiving fiber) so as to achieve a first level of optical signal intensity stabilization at signal detector of the receiver modules.

The object of the Speckle-to-Fiber/Detector Mechanism in the signal reception channels of system is to lock a maximum intensity speckle in the received laser beam onto the receiving fiber, so as achieve a second level of optical signal intensity stabilization at the signal detector in the receiver modules of the system.

FIG. 9E shows a tripod-mounted FS-OA laser communication terminal constructed in accordance with the architecture of FIG. 9B, wherein the transceiver telescope has a 20 cm aperture and is mounted beneath the optical bench of the terminal, and a high-spatial resolution AO (WFS/WFC) compensation subsystem is provided for communication link characterization/compensation, and a real-time subsystem is provided for stabilizing the intensity of laser beam carrier signals detected at the transceiver of the system.

FIG. 9F shows a compact FS-OA laser communication terminal also constructed in accordance with the architecture of FIG. 9B, comprising: a transceiver telescope having a 15 cm aperture and mounted beneath the optical bench of the terminal; a high-spatial resolution AO (WFS/WFC) compensation subsystem with communication link characterization/compensation; a real-time subsystem having a Laser Beam Speckle Tracking Mechanism (including a fast steering mirror FSM)); and a Speckle-to-Fiber/Detector Locking Mechanism (including a deformable mirror as a spatial phase modulation panel) for stabilizing the intensity of laser beam carrier signals detected at the transceiver of the system, in response to atmospheric turbulence.

Free-Space Adaptive Optical (FS-OA) Laser Communication System Supporting Optically-Combined Signal Transmission and Reception Channels and Employing Speckle-To-Receiver-Aperture Locking Mechanism Requiring Receiver-to-Transmitter “Figure of Merit” Communication Feedback, and Laser Beam Speckle Tracking Mechanisms and Speckle-to-Fiber/Detector Locking Mechanisms Along the Signal Reception Channels Thereof for Automatically Stabilizing Variations in the Intensity of Received Laser Beam Carrier Signals Caused by Atmospheric Turbulence Along Said Signal Channels.

In FIGS. 10A and 10B, there is shown a fifth illustrative embodiment of a free-space optical (FSO) laser communication system 40 in accordance with the principles of the present invention, wherein first and second communication terminals 41A and 41B are in communication by way of a super broad-band FSO laser beam communication link 42 having optically-separated signal transmission and reception channels, and wherein the Speckle-To-Receiver-Aperture Locking Mechanism (requiring receiver-to-transmitter “figure of merit” communication feedback) of the present invention is employed along the transmission channels of each communication terminal, and the Laser Beam Speckle Tracking Mechanism and Speckle-to-Fiber/Detector Mechanism of the present invention are employed along the signal reception channels of each communication terminal, for automatically stabilizing variations in the detected intensity of received laser beam carrier signals caused by atmospheric turbulence along said signal channels. As shown in FIG. 10B, each communication terminal has a transmitter (module) 41A and a receiver (module) 41B, each of which will be described in detail below.

The transmitter in each terminal includes a telescopic transmitting aperture, and an optical train embedded with the components of the Speckle-to-Receiver-Aperture Locking Mechanism realized by a fast steering mirror and a quad-cell detector of the kind shown in FIG. 3F, supporting optics, and a processor for carrying out one of the tracking algorithms shown in FIGS. 4A through 4C using the received communication link power signal (sent via the laser communication link), as shown in FIG. 10B and described in detail hereinabove with reference to the description of the corresponding Speckle-to-Receiver-Aperture Locking Method.

The receiver in each terminal includes a telescopic receiving aperture, and an optical train embedded with the components of the two fade mitigation mechanisms of the present invention, namely: the Laser Beam Speckle Tracking Mechanism hereof realized using a FSM, a quad-cell detector shown in FIG. 3F, supporting optics, and a processor carrying out the tracking algorithms shown in FIGS. 4D through 4F; and the Speckle-To-Fiber/Detector Locking Mechanism hereof realized using spatial phase or intensity phase modulator shown in FIG. 3E, a receiving fiber, a single cell detector, and the processor carrying out (i) one of the spatial phase modulation (SPM) control signal generation algorithms shown in FIGS. 5A through 5D when using spatial phase modulation (SPM) techniques, or (ii) one of the spatial intensity modulation (SIM) control signal generation algorithms shown in FIGS. 6A through 6D when using spatial phase modulation (SPM) techniques.

The object of the Laser Beam Speckle Tracking Mechanism employed in both the signal reception channels of system is to automatically track or follow a maximum intensity laser beam speckle and move away from low intensity (i.e. black) laser beam speckles (that might fall onto the receiving fiber) so as to achieve a first level of optical signal intensity stabilization at signal detector of the receiver modules.

The object of the Speckle-to-Fiber/Detector Mechanism in the signal reception channels of system is to lock a maximum intensity speckle in the received laser beam onto the receiving fiber, so as achieve a second level of optical signal intensity stabilization at the signal detector in the receiver modules of the system.

This FSO laser communication system shown in FIGS. 10A and 10B is similar to the laser communication system shown in FIGS. 3A and 3B, with the addition of Speckle-To-Receiver-Aperture Locking Mechanism. As described hereinabove, this is essentially the same basic technique as the Speckle-To-Fiber/Detector Locking Method of the present invention. In this case, however, a spatial phase or intensity modulator is used in the transmitter to optimize the laser beam signal passed into the receiver aperture on the other side of the communication link. To achieve this, it is necessary for the receiver to send the communication link power signal that it is using to perform the Speckle-to-Fiber/Detector Locking Method back to the transmitter via the optical or other communications link. The transmitter uses this signal as the figure of merit in the algorithm that drives its spatial phase or intensity modulator. This control loop acts to maintain (i.e. lock) a bright speckle in the transmitted laser beam onto the receiver aperture, while the receiver loop maintains a bright speckle in the received laser beam in its focal plane onto the receive fiber (or detector).

It is understood that the image capture and processing technology employed in bar code symbol reading systems of the illustrative embodiments may be modified in a variety of ways which will become readily apparent to those skilled in the art of having the benefit of the novel teachings disclosed herein. All such modifications and variations of the illustrative embodiments thereof shall be deemed to be within the scope and spirit of the present invention as defined by the Claims to Invention appended hereto. 

1. In an optical system having-a wave front sensing (WFS) compensation subsystem employing a Hartmann wavefront sensor, in which a stored reference wavefront is represented as a set of reference spot positions on said Hartmann wavefront sensor and a wavefront correction element is employed to correct the wavefront of an optical signal received by a receiver of the optical system or transmitted from a transmitter of said optical system relative to a second optical signal received from the receiver, a wavefront sensing (WFS) control method for controlling the operation of said WFS compensation subsystem comprising the steps of: setting each said reference spot position on said Hartmann wavefront sensor; setting a subaperture region about each said set reference spot position; producing and capturing a plurality of Hartmann spots within said Hartmann wavefront sensor; estimating a distance that each Hartmann spot is from the reference spot of a corresponding subaperture region; generating drive signals for controlling said wavefront correction element, based on the estimated distances; calculating a performance metric that measures a performance of said optical system based on measurements associated with the optical signal; updating said reference spot positions of said stored reference wavefront; and repeating the WFS control method to optimize the performance metric.
 2. The WFS control method of claim 1, which further comprises automatically resetting the subaperture region about each said reference spot position on said Hartmann wavefront sensor when the estimated distances exceed a predetermined threshold.
 3. The WFS control method of claim 1, wherein said optical system is a free-space optical (FSO) laser communication system having a receiving aperture, wherein said optical signal is a laser beam, wherein said FSO laser communication system receives said laser beam using an optical fiber and a signal detector, and wherein said performance metric is the signal power of the detected laser beam.
 4. The WFS control method of claim 1, wherein said optical system comprises a laser beam delivery system having a transmission aperture, wherein said optical signal is a laser beam, wherein said laser beam delivery system transmits said laser beam towards a target, and wherein said performance metric is the signal power of the laser beam delivered to said target.
 5. The WFS control method of claim 1, wherein said optical system comprises a telescopic imaging system having a receiving aperture, wherein said optical signal is a light beam, wherein said telescopic imaging system receives said light beam, and wherein said performance metric is the sharpness of an image formed by said telescopic imaging system.
 6. The WFS control method of claim 1, wherein said wavefront correcting element is a deformable mirror having a plurality of wavefront correcting elements controlled by said generated drive signals.
 7. The WFS control method of claim 1, wherein the steps of producing and capturing the plurality of Hartmann spots, estimating the distance that each Hartmann spot is from the reference spot, generating the drive signals, and calculating the performance metric are performed continuously in a closed-loop manner at a first rate, and wherein the step of updating said reference spot positions is performed in a closed-loop manner at a second rate which is slower than said first rate.
 8. WFS control method of claim 2, wherein the steps of producing and capturing the plurality of Hartmann spots, estimating the distance that each Hartmann spot is from the reference spot, generating the drive signals, and calculating the performance metric are performed continuously in a closed-loop manner at a first rate, and wherein the step of automatically resetting the subaperture region about each said reference spot position is performed in a closed-loop manner at a second rate which is slower than said first rate.
 9. The WFS control method of claim 1, wherein said optical system is a moving platform.
 10. The WFS control method of claim 2, wherein the Hartmann spot positions are averaged over several frames of Hartmann spot data associated with the produced and collected plurality of Hartmann spots so as to determine which Hartmann spots are moving outside of their corresponding subaperture regions; and wherein upon the estimated distances exceeding said predetermined threshold, the subaperture regions about said reference spot positions are automatically reset on said Hartmann wavefront sensor, so as to bring the Hartmann spots closer to the center of their corresponding subaperture regions.
 11. A wavefront sensing (WFS) compensated optical system comprising: a wavefront sensing (WFS) compensation subsystem employing a Hartmann wavefront sensor, in which a stored reference wavefront is represented as a set of reference spot positions on said Hartmann wavefront sensor and a wavefront correction element is employed to correct the wavefront of an optical signal received by a receiver of the optical system or transmitted from a transmitter of said optical system relative to a second optical signal received from the receiver; and a wavefront sensing (WFS) control process operating within said optical system for controlling the operation of said WFS compensation subsystem; wherein said WFS control process comprises the steps of: setting each said reference spot position on said Hartmann wavefront sensor; setting a subaperture region about each said set reference spot position; producing and capturing a plurality of Hartmann spots within said Hartmann wavefront sensor; estimating a distance that each Hartmann spot is from the reference spot of a corresponding subaperture region; generating drive signals for controlling said wavefront correction element, based on the estimated distances; calculating a performance metric that measures a performance of said optical system based on measurements associated with the optical signal; and automatically resetting said reference spot positions of said stored reference wavefront and repeating said WFS control process so as to optimize said performance metric of said optical system.
 12. The WFS-compensated optical system of claim 11, wherein said WFS-control process further comprises automatically resetting the subaperture region about each said reference spot position on said Hartmann wavefront sensor when the estimated distances exceed a predetermined threshold.
 13. The WFS-compensated optical system of claim 11, wherein said optical system is a free-space optical (FSO) laser communication system having a receiving aperture, wherein said optical signal is a laser beam, wherein said FSO laser communication system receives said laser beam using an optical fiber and a signal detector, and wherein said performance metric is the signal power of the detected laser beam.
 14. The WFS-compensated optical system of claim 11, wherein said WFS-compensated optical system comprises a laser beam delivery system having a transmission aperture, wherein said optical signal is a laser beam, wherein said laser beam delivery system transmits said laser beam towards a target, and wherein said performance metric is the signal power of the laser beam delivered to said target.
 15. WFS-compensated optical system of claim 11, wherein said WFS-compensated optical system comprises a telescopic imaging system having a receiving aperture, wherein said optical signal is a light beam, wherein said telescopic imaging system receives said light beam, and wherein said performance metric is the sharpness of an image formed by said telescopic imaging system.
 16. The WFS-compensated optical system of claim 11, wherein said wavefront correcting element is a deformable mirror having a plurality of wavefront correcting elements controlled by said generated drive signals of said WFS control process.
 17. The WFS-compensated optical system of claim 11, wherein the steps of producing and capturing the plurality of Hartmann spots, estimating the distance that each Hartmann spot is from the reference spot, generating the drive signals, and calculating the performance metric in said WSF control process are performed continuously in a closed-loop manner at a first rate, and wherein the step of automatically resetting said reference spot positions is performed in a closed-loop manner at a second rate which is slower than said first rate.
 18. The WFS-compensated optical system of claim 12, wherein the steps of producing and capturing the plurality of Hartmann spots, estimating the distance that each Hartmann spot is from the reference spot, generating the drive signals, and calculating the performance metric in said WFS control process are performed continuously in a closed-loop manner at a first rate, and wherein the step of automatically resetting the subaperture region about each said reference spot position is performed in a closed-loop manner at a second rate which is slower than said first rate.
 19. The WFS-compensated optical system of claim 11, wherein said WFS-compensated optical system is a moving platform.
 20. The WFS-compensated optical system of claim 12, wherein the Hartmann spot positions are averaged over several frames of Hartmann spot data associated with the produced and collected plurality of Hartmann spots so as to determine which Hartmann spots are moving outside of their corresponding subaperture regions; and wherein upon the estimated distances exceeding said predetermined threshold, the subaperture regions about said reference spot positions are automatically reset on said Hartmann wavefront sensor, so as to bring the Hartmann spots closer to the center of their corresponding subaperture regions. 